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ABSTRACT 


The  effect  of  ferromagnetism  on  i ntergranu 1 ar 
segregation  and  solid  solubility  in  Fe  was  theoretically 
investigated  and  the  solid  solubility  of  P  in  alpha  Fe  was 
experimentally  investigated. 

For  the  solubility  measurements,  a  special  apparatus 
was  constructed  to  produce  high  purity  iron  phosphide  from 
elemental  iron  and  phosphorus.  The  iron  phosphide  and  pure 
iron  were  vacuum  arc-melted  to  produce  high  purity  alpha 
Fe-P  alloys.  These  alloys  were  aged  at  temperatures  between 
800  C  and  500  C  and  EPMA  techniques  were  used  to  determine 
the  equilibrium  solubility  of  P.  The  results  agreed  with 
recently  published  data  by  Nishizawa  et  al.  Nishizawa's 
explanation  of  anomalous  solid  solubility  in  alpha  iron 
binaries  below  the  Curie  temperature  by  considering 
non-linear  magnetic  components  of  thermodynamic  functions 
was  confirmed. 

In  the  theoretical  part  of  the  work  empirical  magnetic 
components  of  thermodynamic  functions  for  alpha  iron  were 
calculated  from  data  in  the  literature.  These  thermodynamic 
functions  were  applied  to  confirm  Nishizawa's  predicted 
behavior  of  solute  solubility  in  an  alpha  iron  binary.  The 
formalism  was  then  extended  to  predict  equilibrium  vacancy 
concentration  in  alpha  iron.  In  addition  Guttmann's  approach 
to  grain  boundary  solute  segregation  for  alpha  iron  binaries 
and  ternaries  was  treated  using  magnetic  components.  It  is 
thereby  predicted  that  most  non-magnetic  solute  in  iron  is 
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expected  to  segregate  to  grain  boundaries  in  larger  amounts 
at  higher  temperatures  than  predicted  by  a  purely 
non-magnetic  analysis.  Unexpected  magnetic  transitions 
within  the  grain  boundary  are  also  predicted  to  occur. 

A  consequence  of  this  treatment  is  that  because  of  its 
ferromagnetic  nature,  alpha  iron  is  more  susceptible  to 
intergranular  segregation  and  thus  intergranular 
embrittlement  than  an  equivalent  non-magnetic  alloy. 
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1.  INTRODUCTION 


Two  processes  which  can  act  simultaneously  in  an  iron 
alloy  are  grain  boundary  solute  segregation  and  second  phase 
precipitation.  Both  can  affect  the  structure  and  certain 
properties  of  the  alloy.  One  such  property  is  intergranular 
cohesion,  a  decrease  in  which  can  result  in  intergranular 
fracture.  In  many  iron  alloys  this  intergranular 
embrittlement  is  associated  with  certain  heat  treatments 
(causing  temper  embrittlement)  and  is  generally  manifest  as 
a  reduction  in  impact  strength  and  fracture  toughness  [1]. 

Segregation  is  defined  as  the  non-uniform  distribution 
of  elements  in  a  material.  Frequently  used  in  reference  to 
castings,  where  segregation  occurs  during  solidification  of 
a  melt,  intergranular  segregation  in  the  present  context 
refers  to  the  solid  state  thermal  redistribution  of  species 
between  the  matrix  and  interfaces.  It  is  also  necessary  to 
distinguish  between  two  types  of  solid  state  segregation 
phenomena.  A  non-equilibrium  type  of  segregation  may  be 
connected  with  the  dragging  of  solute  species  by  a  flux  of 
vacancies  toward  the  interface  [2]  during  cooling  from  a 
high  temperature  or  during  irradiation  [3].  Another 
non-equilibrium  type  of  segregation  involves  the  rejection 
of  solute  during  precipitate  growth. 

However,  the  omnipresent  form  of  i ntergranu 1 ar 
segregation  is  a  reversible  equilibrium  segregation  [4].  For 
a  polycrystalline  alloy,  held  i sotherma 1 ly ,  there  will  be  an 
enrichment  of  surface  active  species  at  interfaces.  The 
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extent  of  this  enrichment  is  determined  only  by  the  system 
parameters  at  equilibrium  and  not  by  their  history. 
Therefore,  as  the  temperature  of  the  alloy  is  changed,  there 
are  compositional  and  perhaps  structural  changes  at  the 
grain  boundaries.  From  this  point  the  word  segregation  will 
be  used  to  denote  this  phenomenon. 

The  precipitation  process  in  a  material  requires  the 
nucleation  and  growth  of  a  second  phase  from  a 
super-saturated  solid  solution.  The  second  phase  may  be 
nucleated  heterogeneously  at  grain  boundaries  and 
interfaces,  dislocations  and  other  structural  defects,  or 
homogeneously  within  the  matrix.  However,  at  equilibrium  for 
a  given  temperature,  there  is  a  limit  of  solubility  for  each 
solute  species  present  in  the  material. 

The  present  project  stems  from  the  postulate  that 
precipitation  and  i ntergr anu 1 ar  segregation  in  iron  must  be 
simultaneous  consequences  of  the  tendency  of  the  iron  to 
reject  foreign  atoms  from  its  matrix.  This  rejection  of 
solute  atoms,  which  is  reflected  in  the  limit  of  solid 
solubility  and  in  grain  boundary  segregation,  is  best 
treated  by  applying  conditions  of  thermodynamic  equilibrium 
to  the  system.  Of  special  interest  in  iron  alloys  are  the  so 
called  'impurity  elements'  ( S , P , As , Sb , Sn  etc.)  which  are 
highly  surface  active  and  are  associated  with  temper 
embrittlement  phenomena. 

It  is  well  Known  that  the  solubility  of  element  X  is 
very  sensitive  to  the  relative  position  of  the  free  energy 
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curve  of  the  precipitate  (containing  species  X)  with  respect 
to  the  free  energy  curve  of  the  matrix.  Guttmann  [5]  has 
recently  shown  that  equilibrium  grain  boundary  segregation 
is  also  very  sensitive  to  the  relative  position  of  the  free 
energy  curve  of  the  segregated  interface  with  respect  to  the 
free  energy  curves  of  the  bulk  phases.  This  correlation 
between  solubility  and  intergranular  segregation  implies 
that  the  energetic  properties  of  a  boundary  are  related  to 
the  energetic  properties  of  a  1 1  the  matrix  phases. 

Therefore,  material  properties  and  external  conditions 
affecting  precipitation  in  a  given  material  should  also 
affect  grain  boundary  segregation  although  the  exact  nature 
of  the  effect  is  not  obvious.  This  is  very  important  in  that 
data  on  precipitation  in  simple  systems  is  relatively  easily 
available  while  good  segregation  data  for  simple  systems  is 
difficult  to  obtain  and  consequently  scarce. 

When  specifically  considering  iron,  the  one  obvious  but 
often  ignored  material  property  is  its  ferromagnetic  nature. 
The  following  3  paragraphs  have  been  compiled  from 
references  [6-8]  as  a  brief  description  of  ferromagnetism  in 
i  ron . 

The  3d  shell  of  the  Fe  atom  is  responsible  for  its 
ferromagnet ic  behaviour.  The  3d  shell  may  have  up  to  10 
electrons,  gradually  filling  from  scandium  to  copper  across 
the  first  row  of  transition  elements.  For  each  electron 
within  the  shell  and  along  some  specified  direction,  there 
is  a  component  of  orbital  angular  momentum  which  must  be  an 


, 
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integral  multiple  of  h/27T.  In  addition,  each  electron  has  an 
intrinsic  angular  momentum,  called  spin,  due  to  a  rotation 
about  its  own  axis;  this  spin  can  be  either  +1/2(h/2^)  or 
-1/2(h/2Tr)  along  the  specified  direction.  The  spin  of  a 
charged  particle  is  associated  with  a  magnetic  moment.  A 
fundamental  principle  of  atomic  structure,  Pauli's  exclusion 
principle,  states  that  no  two  electrons  in  a  given  shell  may 
have  identical  components  of  orbital  angular  momentum  with 
identical  components  of  spin.  However,  the  manner  in  which 
electrons  are  added  to  a  shell  as  the  atomic  number 
increases  is  important.  Hund' s  rule  states  that  within  a 
given  shell,  the  maximum  number  of  electrons  have  spins 
pointing  in  the  same  direction,  consistent  with  Pauli's 
exclusion  principle.  There  is  a  coupling  of  exchange  energy 
between  electrons  with  parallel  spins;  this  exchange  energy 
is  negative  for  electrons  in  the  same  atom  and  leads  to  a 
minimum  energy  when  the  maximum  number  of  spins  are 
paral lei . 

When  a  number  of  atoms  are  associated,  as  in  a  crystal, 
this  internal  interaction  is  called  the  exchange  field, 
molecular  field,  or  Weiss  field  (after  Pierre  Weiss).  The 
exchange  field  gives  an  appropriate  representation  of  the 
quantum  mechanical  exchange  interaction.  It  has  been  shown 
that  the  energy  of  interaction  w^ 
spins  s i ,  s^  contains  a  term 


of  atoms  1,  2  bearing 
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wi 2 =  - j  srs2 

(Hei senberg  model ) 

(1.1) 

where  j  is  the  exchange  integral  and  is  related  to  the 
overlap  of  the  charge  distribution  of  the  atoms. 

The  outer  shells  of  an  isolated  Fe  atom  are  3d64s2  and 
the  3d  shell  contains  4  coupled  electrons  (antiparallel 
electrons  cancel  each  other).  In  general,  the  parallel  spin 
arrangement  will  be  energetically  favoured  over  the 
antiparallel  one.  This  condition  exists  in  the  electronic 
(conduction)  band  of  iron,  cobalt  and  nickel.  The  density  of 
levels  within  the  band  is  so  high  that  the  negligible  gain 
in  kinetic  energy  required  for  an  electron  to  move  from  a 
filled  level  to  an  unfilled  one  is  more  than  offset  by  the 
lowering  of  electrostatic  energy.  This  is  responsible  for 
the  alignment  of  electron-spin  magnetic  moments  that 
constitute  the  spontaneous  magnetization  existing  within 
domains  of  ferromagnetic  materials.  Below  a  critical 
temperature  (material  dependent)  called  the  Curie 
temperature,  there  is  a  ferromagnet i c  coupling  between 
electrons  with  parallel  spins,  but  above  this  temperature, 
the  coupling  is  destroyed  by  thermal  energy  (paramagnetism). 

The  above  description  is  included  for  completeness  but 
it  is  not  the  intent  of  the  present  work  to  delve  further 
into  the  quantum  mechanical  aspects  of  ferromagnetism. 
Rather,  an  equilibrium,  macroscopic  thermodynamic  approach 
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will  be  considered.  Empirical  thermodynamic  functions 
obtained  from  solid  solubility  data  in  alpha  iron  will 
used  to  predict  grain  boundary  segregation  behavior  of 
selected  solute  species  in  alpha  iron. 


be 


2.  BACKGROUND 

The  purpose  of  this  chapter  is  to  provide  a  brief 
historical  background  to  the  present  work.  The  development 
of  grain  boundary  segregation  theory  is  presented  with  some 
reference  to  the  segregation  of  phosphorus  in  alpha  iron. 
Excellent  reviews  of  grain  boundary  segregation  have  been 
published  by  Inman  and  Tipler  [9]  and  by  Hondros  and  Seah 
[4].  As  the  experimental  determination  of  phosphorus 
solubility  in  alpha  iron  was  part  of  this  project,  a  brief 
review  of  previous  investigations  into  P  solubility  is  also 
presented . 


2.1  The  Development  of  Segregation  Theory 

The  driving  force  for  interfacial  segregation  must  be 
the  lowering  of  interfacial  energy.  The  theory  applies  to 
all  interfaces  but  the  major  interest,  at  present,  is  in 
grain  boundaries.  Grain  boundaries,  which  form  due  to  the 
geometrical  discontinuity  between  two  adjacent  lattices,  are 
distorted  and  more  or  less  disordered  areas  where  the 
orientation  and  maqni tude  of  atomic  bonds  are  different  from 
those  in  the  bulk  crystal.  It  can  be  considered  to  be  a 
discontinuity  large  in  two-dimensional  extent  but  only  a  few 
atomic  diameters  thick.  Within  this  volume,  it  can  be 
treated  as  distorted  lattice  containing  compressed  and 
expanded  lattice  sites. 
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About  one  hundred  years  ago  Gibbs  [10]  developed  a 
thermodynamic  treatment  of  surfaces  that  included  such 
considerations  as  adsorption,  heat  and  mechanical  surface 
effects,  and  grain  boundaries.  This  is  the  framework  for  all 
modern  thermodynamic  theory.  For  a  surface  in  a  binary 
system  with  components  1  and  2,  Gibbs  derived  the  following: 


( dy/ dy i )T^p 


[r2  -  x^] 


x 


(2.1) 


where  y  is  the  surface  energy  of  the  interface, 
jj-L  is  the  chemical  potential  of  1, 

X.  is  the  bulk  atomic  fraction  of  i,  and 
r •  is  the  excess  surface  concentration  of  i. 

This  familiar  Gibbs  adsorption  isotherm  is  the  starting 
point  of  segregation  theory. 

The  earliest  segregation  calculations  were  made  in  the 
development  of  a  theory  of  yield-point  behavior  in  metals  by 
Cottrell  [11].  He  considered  that  the  binding  energy,  E, 
between  a  solute  atom  and  a  dislocation  could  be  evaluated 
by  considering  the  release  of  strain  energy  on  segregation 
to  the  dislocation1,  such  that  in  a  dilute  case, 


C  =  CQexp(E/RT) 


(2.2) 


where  C  is  the  local  concentration  of  the  solute  and 
Co  is  the  average  bulk  concentration  of  solute. 


’For  example,  E  for  nitrogen  in  alpha  iron  was  estimated  at 
50  KJ/mole. 
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In  the  mid  1950's  McLean  [12]  considered  equilibrium 
segregation  in  terms  of  the  lattice  distortion  energy  around 
solute  atoms  in  the  crystal  bulk.  A  solute  atom  larger  than 
the  site  it  occupies  in  the  solvent  lattice  should  fit 
easily  in  an  expanded  boundary  site  causing  less  lattice 
distortion.  Likewise,  a  smaller  solute  atom  should  fit  more 
easily  into  a  compressed  site.  Generally,  large  energies  are 
required  to  alter  the  volume  occupied  by  atoms  in  a  metal. 
For  example,  an  elastic  10%  expansion  in  the  alpha  iron 
lattice  requires  approximately  62.7  KJ/mole.  The  solution  of 
atoms  different  in  size  from  the  solvent  requires  strain 
energy  which  can  be  relieved  by  segregation  of  the  solute  to 
distorted  sites.  In  this  way,  a  large  amount  of  the 
distortion  energy  of  the  entire  system  can  be  released  by 
segregat ion . 

McLean  produced  the  useful  expression 

X41  =  exp  (  E/RT) 

1-X4  1-Xa 

(2.3) 

where  X4  =  the  boundary  solute  content 
Xa  =  the  bulk  solute  content 

E  =  the  binding  energy  of  the  solute.  This  model  shows 
that  segregation  rises  as  the  bulk  solute  content  rises  or 
as  the  temperature  falls.  Segregation  will  occur  in  nearly 
all  solid  solutions.  However,  it  would  be  noticeable  only 
when  the  temperature , t ime  and  diffusion  coefficient  are  such 
that  a  large  grain  boundary  concentration  can  build  up.  To 
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grain  boundary  segregation  it  is  necessary  to  quench  from  a 
temperature  high  enough  to  disperse  solute  to  one  low  enough 
to  prevent  significant  diffusion. 

McLean's  model  has  the  same  behavior  as  described  by 
the  simple  Langmuir  theory  of  gas  adsorption  on  a  solid 
surface.  That  is,  there  is  a  single  layer,  a  single 
adsorbate  and  no  site-to-site  interaction.  There  are  very 
few  grain  boundary  results  for  iron  binaries  that  are 
sufficiently  complete  to  be  properly  interpreted  in  terms  of 
McLean- Langmui r  type  theories  [13,14]. 

Hondros  and  Seah  [13]  have  invoked  other  solid  state 
analogues  of  theories  developed  for  gas  adsorption  on  solid 
surfaces,  adopting  the  concepts  and  analytical  treatments 
commonly  used  in  classical  free  surface  adsorption.  One  of 
these  is  the  BET  theory  of  multilayer  gas  adsorption  [15], 
given  by 


X*  X“  1  K-l  Xa 

x1*1  x“-xa  K  K  x“ 

(2.4) 


where  K 
E 


y  $ 
*1 

Xa 


exp( E/RT ) 

the  difference  between  the  energy 
of  adsorption  of  the  first  layer  and 
the  free  energy  of  condensation  of 
successive  layers 
boundary  solute  content 
boundary  solute  content  of  a  layer 
bulk  solute  content 


1 1 


X1^  =  maximum  bulk  solute  content 
In  this  model  it  was  assumed  that  atoms  adsorbing  in  the 
first  monatomic  layer  have  a  different  free  energy  of 
adsorption  than  those  in  successive  layers.  The  adsorbate 
atoms  have  no  lateral  interaction  in  the  layer  and  each 
layer  does  not  have  to  be  complete  before  the  next  is 
started.  A  truncated  BET  model2  is  very  similar  to  the 
McLean  model . 

The  above  theories  neglect  site-to-site  interactions 
and  chemical  interactions.  Site-to-site  interactions  are 
mutual  interactions  between  the  same  solute  in  neighboring 
boundary  sites  and  chemical  interactions  are  element-element 
interactions  found  in  systems  with  multicomponent 
segregation.  In  the  Fowler  model  [16]  there  is  a  single 
layer,  a  single  segregating  solute  species,  and  site-to-site 
interaction.  For  the  gas  adsorption  models,  the  assumption 
that  a  filled  site  does  not  affect  the  probability  of  one  of 
z  neighbouring  sites  being  filled  can  be  removed  by  allowing 
an  interaction  energy  r  between  pairs  of  nearest  neighbours. 
The  solid  state  analogue  is 

_X*  =  Xaexp({E-rz[X^/X^]}/RT) 

1'X<(>  (2.5) 

If  r  is  positive  the  atoms  are  mutually  repulsive  and  if 
negative  vice  versa.  With  r  negative  the  temperature 

2The  boundary  conditions  are  the  same  as  in  the  McLean 
mode  1 . 
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dependence  of  segregation  deviates  from  that  given  by  McLean 
by  showing  a  sharper  rise  in  segregation  as  temperature 
falls. 

In  the  recent  work,  Guttmann  [17]  has  developed  a 
thermodynamic  treatment  of  mul t i -component  segregation.  His 
model  is  used  in  the  present  work  and  it  will  therefore  be 
discussed  fully  in  chapter  4. 

The  particular  system  of  major  interest  here  is  alpha 
Fe-P.  There  has  been  a  great  deal  of  experimental  work 
concerning  P  segregation  in  steels3  as  P  is  a  potent 
i ntergranul ar  embrittling  element,  but  little  on  pure  alpha 
Fe-P  binaries.  The  reason  is  that  it  is  difficult  to  temper 
embrittle  pure  iron  binaries  sufficiently  so  that 
i ntergranu 1 ar  fracture  may  occur  and  the  surface  may  be 
studied.  In  the  late  1950's  Inman  and  Tipler  [18]  carried 
out  a  study  of  P  segregation  in  Fe  by  using  radioactive  P. 
The  occurrence  of  segregation  was  conclusively  shown  but  the 
resolution  was  poor  since  a  sectioning  technique  was  used. 
More  recent  investigations  [14,19]  have  used  Auger  Electron 
Spectrometer  (AES)  techniques.  However,  the  results  are 
scattered  with  widely  varying  saturation  levels  and  do  not 
seem  to  fit  a  simple  McLean- type  model. 


3 For  example,  see  the  work  of  McMahon  at  the  University  of 
Pennsy 1 vani a . 
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2.2  The  Solubility  of  Phosphorus  in  alpha  Iron 

Alpha  iron  forms  a  substitutional  solid  solution  with 
phosphorus  to  a  maximum  solubility  of  2.55  wt%  at  -1050  C. 
Phase  diagrams  and  a  compilation  of  early  works  are  given  by 
Hansen  [20]  and  Vol  [21];  these  demonstrate  an  interest  in 
the  Fe-P  system  since  before  the  turn  of  the  century. 

The  solid  solubility  of  P  in  alpha  Fe  decreases  with 
decreasing  temperature.  Before  the  advent  of  quantitative 
Electron  Probe  Microanalysis  (EPMA),  the  most  accurate 
solubility  determinations  were  by  meta 1 lographic  analyses4 
[22,23].  Temperatures  from  1000  C  to  600  C  were  investigated 
but  the  most  reliable  values  lie  above  the  Curie  temperature 
of  i ron  (  1  043  K )  . 

In  the  1960's  Hornbogen  [24-26]  investigated  the 

factors  that  govern  precipitation  from  alpha  Fe 

substitutional  solid  solutions.  Phosphorus  was  chosen 

because  of  its  small  atomic  size  as  compared  to  iron  and  the 

lack  of  chemical  and  crysta 1 lographi c  similarity  between  the 

two  elements.  He  reported  that  high  P  alloys  tended  to  show 

extensive  i ntergr anu 1 ar  brittleness.  The  phosphorus 

precipitated  from  alpha  Fe  as  Fe^P  when  equilibrium  was 

reached.  This  phosphide  which  has  a  tetragonal  body-centered 

unit  cell  containing  32  atoms,  forms  precipitates  parallel 

to  a  high  index  plane  of  alpha  Fe  (1,4,12).  Hornbogen  also 

found  that  phosphides  nucleate  from  alpha  iron  in  the 

4Resistivity  measurements  and  lattice  parameter 
determinations  by  X-ray  techniques  were  not  very  successful 
[23]  . 
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sequence:  at  high  angle  boundaries,  at  low  angle  boundaries 

and  dislocations,  and  in  the  matrix. 

More  recent  determinations  of  P  solubility  [27-29]  have 

involved  the  use  of  EPMA  and  are  shown  in  Figure  2.1. 

Hofmann  et  al.  [27]  were  among  the  first  to  notice  an 

anomalous  decrease  in  the  solubility  of  P  and  in  its 

diffusion  coefficient  in  iron  below  the  Curie  temperature. 

More  recently,  Doan  et  al.  [28]  in  investigating  the  Fe-Ni-P 

system  noticed  a  much  more  pronounced  anomalous  decrease  in 

P  solubility.  While  the  present  project  was  in  progress, 

Nishizawa  et  al.  [29]  published  an  even  more  striking 

decrease  in  P  solubility  below  the  Curie  temperature,  and 

related  it  to  the  ferromagnetism  of  alpha  iron  (see 

Chapter  4).  It  is  very  interesting  that,  as  EPMA  techniques 
% 

have  improved  and  more  time  allowed  for  alloys  at 
temperature  to  reach  equilibrium,  the  experimentally 
measured  solubility  of  P  in  alpha  iron  has  been  continously 
decreas i ng . 


Solubility  of  P  in  bcc  Ft  /  at% 
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Figure  2.1  The  solubility  of  phosphorus  in  alpha  iron 


3.  EXPERIMENTAL  PROCEDURES 

Although  this  thesis  is  primarily  theoretical  some  initial 
experimental  work  was  performed.  This  experimental  work  was 
quite  difficult  and  so  the  procedures  used  are  described  in 
extensive  detail  in  this  section. 


3.1  The  Production  of  'Pure'  alpha  Fe-P  Alloys 

The  production  of  'pure'  alloys  containing  only  iron 
and  phosphorus  was  necessary  in  order  to  determine  the 
equilibrium  limit  of  solid  solubility  of  phosphorus  in  alpha 
iron.  The  presence  of  additional  elements  could  affect  the 
behavior  of  the  alloys.  For  example,  at  higher  temperature, 
small  quantities  of  a  similar  atom  such  as  Ni  are  tolerable 
for  experimental  purposes  whereas  small  quantities  of  Mo  are 
not.  This  is  illustrated  by  some  of  the  results  obtained 
(see  e.g.  section  4.11).  The  first  step  in  the  production  of 
the  'pure'  alpha  Fe-P  alloys  must  be  the  procurement  of  the 
proper  raw  materials. 

3.1.1  Sources  of  Iron  and  Phosphorus 

An  initial  decision  that  was  faced  was  whether  to  use 
phosphorus  in  the  elemental  form  or  as  a  compound.  In  the 
elemental  form,  phosphorus  exists  in  a  number  of  allotropes 
including  white  (yellow),  red,  and  black  (violet).  White 
phosphorus  burns  spontaneously  in  air  producing  a  pentoxide 
which  is  very  poisonous  (50  mg  fatal  dose).  However,  when 
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heated  in  its  own  vapour  to  250  C,  white  phosphorus  is 
converted  to  the  red  variety.  This  form  is  fairly  stable, 
does  not  ignite  in  air  but  sublimes  with  a  vapour  pressure 
of  1  atm  at  417  C.  For  the  present  project,  red  phosphorus 
(99.9+  wt%)  was  obtained  from  Fisher  Scientific  Co.  This 
material  was  Kept  sealed  under  argon  gas  as  much  as  possible 
to  minimize  oxygen  contamination. 

The  only  phosphorus-containing  compounds  suitable  for 
production  of  alloys  for  this  project  are  iron  phosphides. 
Ferrophosphorus  is  used  as  an  alloying  addition  by  the  steel 
industry.  This  material  has  a  high  phosphorus  content 
(—25  wt%);  however,  it  also  contains  varying  and 
unacceptable  amounts  of  carbon,  silicon,  manganese  and 
chromium.  Commercially,  there  are  very  few  sources  of 
phosphides  available.  ICN  Pharmaceuticals,  Inc.  sells  a 
ferric  phosphide,  FeP,  but  only  guarantees  purity  to  greater 
than  95  wt%.  A  quantity  of  this  commercial  FeP  was  obtained; 
however,  upon  obtaining  a  qualitative  X-ray  spectrum  of  this 
material  using  an  Energy  Dispersive  Spectrometer  (EDS), 
peaks  denoting  significant  quantities  of  Cr  and  Mn  were 
clearly  visible.  Therefore,  this  'pure'  commercial  iron 
phosphide  was  not  used.  Instead  'pure'  iron  phosphide  was 
produced  from  'pure'  iron  and  'pure'  phosphorus. 

Several  commercial  sources  of  'pure'  iron  are 
available.  Three  materials  were  considered  as  the  possible 
iron  base  for  the  alpha  Fe-P  alloys: 

1.  High  Purity  Armco  iron  in  the  form  of  0.017  in.  thick 
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sheet  was  tried,  but  rejected  due  to  the  impurities 
present.  For  example,  Plate  3.1  shows  an  X-ray  spectrum 
from  one  of  many  inclusions. 

2.  Analysed  reagent  grade  iron  wire  was  obtained  from 
J.  T.  Baker  Chemical  Co.  This  material  was  somewhat 
better  than  Armco  iron  and  was  used  for  preliminary 
tests  and  initial  production  of  iron  phosphides. 

3.  Specpure  iron  0.25  in.  diameter  rod  ( spectrographica 1 ly 
standardized)  was  obtained  from  Johnson  Matthey 
Chemicals  Limited.  This  material  is  extremely  pure  and 
expensive  and  therefore  was  used  only  for  the  final 
samples . 

Table  3.1  contains  the  compositions  of  the  above  iron 
mater i a  1 s . 

3.1.2  Production  of  an  Iron  Phosphide 

There  are  several  possible  techniques  for  producing 
iron  phosphide  of  a  desired  composition  from  iron  and 
phosphorus,  including  the  following: 

1.  A  special  vacuum  furnace  may  be  constructed  such  that 
above  a  particular  temperature  pure  iron  is  surrounded 
by  and  thus  reacts  with  phosphorus  vapour.  The  inner 
lining  of  the  furnace  must  resist  attack  from  the 
phosphorus  vapour.  There  is  normally  a  large  loss  of  P 
with  this  met  hod . 

2.  An  ordinary  vacuum  furnace  may  be  used  if  the  amount  of 

P  is  small  compared  to  the  iron.  This  allows  the  P  to  be 
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Table  3.1  Compositions  of  'pure' 

irons  used 

( in  wt%) 

Element 

1 .  Armco 

2.  Wire 

3 .  Specpu 

Ag 

1  ppm 

Ca 

1  ppm 

Mg 

1  ppm 

C 

0.007 

- 

Mn 

0.010 

1  ppm 

P 

0.020 

0.010 

- 

S 

0.002 

0.004 

- 

Cu 

0.130 

<  1  ppm 

Ni 

0.042 

- 

Cr 

0.003 

<  1  ppm 

Cb 

0.005 

- 

Mo 

0.006 

- 

Si 

0.001 

0.005 

<  1  ppm 

Co 

0.009 

- 

Pb 

0.001 

- 

Sn 

0.010 

- 

A 1 

0.003 

- 

Fe 

ba 1 ance 

99.9+ 

ba 1 ance 

A  blank  signifies  the  element  was  analysed  for; 
whereas  indicates  a  concentration  below  the 
limit  of  detectability. 
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Plate  3.1  An  X-ray  spectrum  from  an  inclusion  in  Armco  iron 
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enclosed  within  the  iron  and  to  be  Kept  away  from  the 
furnace  lining. 

In  the  present  work  a  master  alloy  (nominally  20  wt%  P) 
was  produced  by  reacting  iron  with  red  phosphorus  in  a 
weight  ratio  of  4  to  1.  This  was  carried  out  in  the  special 
apparatus  shown  in  Figure  3.1.  A  pyrex  glass  tube,  see 
Plate  3.2,  was  fitted  with  a  fretted  quartz  glass  disc. 

Iron,  in  the  form  of  powder,  turnings,  or  sectioned  sheet 
and  rod  was  placed  inside  the  tube  above  the  disc.  When  the 
tube  was  in  a  vertical  position,  the  disc  Kept  solid 
material  at  a  specific  position  (reaction  chamber)  but 
allowed  passage  of  gases.  The  tube  was  placed  within  a 
vertical  furnace,  see  Plate  3.3.  The  bottom  of  the  tube 
(below  the  furnace)  was  connected  by  a  2-way  valve  to  a 
cylinder  of  ultra-high  purity  hydrogen  gas  and  a 
rotary-diffusion  vacuum  pump  combination. 

The  top  of  the  tube  (above  the  furnace)  was  also 
connected  by  a  2-way  valve  to  a  bubbler  and  to  an  upper 
container  of  red  phosphorus  which  was  directly  connected  to 
an  argon  cylinder.  The  pyrex  tubes  had  a  reduced  internal 
diameter  on  either  side  of  the  reaction  chamber  to  allow  for 
easy  sealing  of  the  tube  with  a  torch. 

The  use  of  the  apparatus  is  summarized  as  follows: 

1.  Iron  was  first  pickled  in  HC1  to  remove  surface 
contamination  and  placed  in  the  reaction  chamber. 
Phosphorus  was  positioned  in  the  upper  chamber  under 
argon.  The  tube  was  then  inserted  in  the  furnace. 
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Figure  3.1  A  schematic  of  the  apparatus  developed  for 
iron  phosphide  production. 
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Plate  3.2  Pyrex  reaction  tube  for  producing  iron  phosphides 


Plate  3.3  Experimental  apparatus  for  producing  iron 
phosphides 
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2.  Hydrogen  gas  was  passed  through  the  tube  to  the  bubbler 
at  approximately  1  bubble  per  second.  The  furnace  was 
then  raised  to  585  C  and  held  for  up  to  24  hours.  This 
procedure  was  to  reduce  iron  oxides,  eliminating  oxygen 
and  other  impurities  [30]. 

3.  When  the  tube  had  cooled  to  room  temperature,  the 
hydrogen  gas  was  shut  off,  the  top  valve  closed  and  the 
tube  evacuated  by  the  pump.  By  opening  the  top  valve  to 
the  phosphorus,  the  powder  was  'sucked'  into  the 
reaction  chamber. 

4.  The  evacuated  reaction  chamber  (tube)  was  then  sealed  by 
torch  at  the  reduced  section. 

5.  The  tube  was  removed  from  the  vertical  furnace,  checked 
for  vacuum  integrity  using  a  Tesla  coil,  and  shaken  to 
mix  the  iron  and  phosphorus  contents. 

6.  The  tube  was  then  placed  completely  within  a  horizontal 
furnace,  slowly  raised  to  585  C  and  held  for  8  hours.5 
During  this  treatment  the  iron  and  phosphorus  reacted  to 
giving  mainly  FeP  with  some  remaining  Fe. 

7.  When  cool,  the  tube  was  carefully  broken  under  an  inert 
atmosphere  and  the  iron  and  phosphide  removed.  Each  tube 
was  thus  used  only  once. 

Although  the  phosphorus  has  a  very  high  vapour  pressure 
and  the  pyrex  tube  is  at  the  point  of  softening  at  the 
temperatures  mentioned,  the  glass  tube  did  not  bulge  or 

5 A 1  though  Fe-P  alloys  have  been  sintered  in  the  past  [31], 
the  present  conditions  were  determined  by  trial. 
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explode  because  the  iron  and  phosphorus  combined  very 
rapidly,  preventing  the  occurrence  of  high  phosphorus  vapor 
pressure . 

Some  of  the  produced  phosphide  was  ground  to  a  powder 
and  analysed  by  X-ray  diffraction  methods.  This  showed  that 
most  of  the  phosphorus  was  present  as  FeP  with  minor  amounts 
of  Fe^P.  It  also  became  apparent  that  unless  the  phosphide 
powder  was  Kept  under  an  inert  atmosphere,  it  oxidized.  An 
oxide  was  particularly  noticeable  when  the  phosphide  powder 
was  vacuum  arc-melted  into  'buttons' .  However,  if  melted 
into  buttons  immediately,  only  stable  phophides  were  present 
and  oxidation  was  not  a  problem. 

3.1.3  Production  of  Final  alpha  Fe-P  Alloys 

All  final  phosphides  (master  alloys)  and  alpha  Fe-P 
alloys  were  vacuum  arc-melted  into  'buttons'  using  a 
Materials  Research  Corporation  V-4  Series  vacuum  arc 
furnace,  Plate  3.4.  A  partial  vacuum  of  argon  was  used.  The 
gas  passed  through  a  dessicator,  over  hot  Cu  turnings 
(600  C),  and  through  a  cold  trap  before  reaching  the  alloys. 
Melting  was  carried  out  on  a  water-cooled  copper  hearth 
using  a  tungsten  rod  electrode.  All  gases  from  the  vacuum 
chamber  were  vented  to  a  fume  hood  to  prevent  accidental 
exposure  to  phosphorus  pentoxides.  Figure  3.2  shows  a 
schematic  of  experimental  procedures  used  in  alpha  Fe-P 
alloy  production.  As  the  figure  shows,  the  Specpure  iron 
rods  were  machined  into  cylinders  2.50  cm.  long,  a  hole 
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Specpure  Iron, 

Va"  Rod  Machined  Into 


Red  Phosphorus 
Power  (  P4 ) 


P4  Sealed  in  Fe  -« 


Plug  with 
5  Degree 
Taper 


Phosphide 

Production 


Iron  Turnings  in 
Reaction  Chamber 


24  hrs  in  H2  Gas 
at  585=0 


Red  Phosphorus 
Power  (  P4 ) 


Method  1  Method  2  Method  3 

(preferred) 

Fe  +  P4 
Vacuum  Sealed 
in  Quartz 


Sintered  in  8  hrs 
at  >585°C 


Fe  +  FeP 


Mixed  and 
Vacuum  Sealed 

Sintered  for  8  hrs 
at  585°C 


Fe  +  FeP 
or 


Vacuum  arc-melted  Stored  Under  Ar 

i 

Fe3P  +  Fe2P  ‘Button’ 


Fe  +  Phosphide  — 

i 

Vacuum  arc-melted 

i 

Final  aFe-P  Alloy 
‘Buttons’ 


Figure  3.2  A  schematic  of  procedures  developed  for  the 
production  of  aFe— P  alloys. 


Plate  3.4  The  vacuum  arc  furnace  used  to  produce  final 
a  1  loys 
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0.318  cm.  i.d.  was  drilled  into  each  cylinder  to  contain 
either  phosphorus  or  phosphide,  and  2.50  cm.  long  plugs  with 
a  5  degree  locking  taper  were  used  to  seal  the  compartment. 
The  cylinders  and  plugs  are  shown  in  Plate  3.5.  These 
configurations  allowed  three  successful  methods  of  final 
alloy  production : 

1.  Red  phosphorus  was  enclosed  within  the  iron  cylinders 
and  directly  vacuum  arc-melted.  This  method  produced 
clean  alloys  but  the  phosphorus  had  to  be  packed  and 
sealed  under  an  inert  atmosphere.  A  limited  amount  of 
powder  could  be  packed  and  at  least  1/3  of  the 
phosphorus  was  lost  on  melting  (larger  cylinders  were 
difficult  to  handle  and  melt).  This  resulted  in  alloys 
containing  less  than  1  wt%  P. 

2.  This  method  was  identical  to  the  above  except  that  the 
sealed  iron  cylinders  were  also  encapsulated  in  quartz 
and  annealed  at  greater  than  600  C  for  24  hours  before 
melting.  This  resulted  in  a  lower  phosphorus  loss. 

3.  The  method  preferred  and  the  one  actually  used  required 
that  the  turnings  from  the  machining  of  the  Specpure 
iron  rods  were  produced  into  phosphide  buttons  as 
described  in  section  3.1.2.  Alloys  of  any  phosphorus 
concentration  could  then  be  made  simply  by  arc-melting  a 
portion  of  the  phosphide  button  with  the  Specpure  iron 
rod  and  atmosphere  was  not  a  problem. 

All  alpha  Fe-P  specimens  produced  were  nominally  1.5  wt%  P. 
The  specimens  were  marked  by  a  letter  and  a  temperature, 
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Plate  3.5  Specpure  iron  cylinders  and  plugs  with  locking 
taper 
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e.g.  T800  C.  T  was  the  designation  used  for  the  final, 
completely  successful  series  using  Specpure  iron. 


3.2  Heat  Treatment  of  alpha  Fe-P  Alloys 

As  previously  stated,  the  solubility  limit  of  P  in  Fe 
above  the  Curie  temperature  (770  C)  has  been  well 
documented.  Therefore,  an  investigation  range  of  800  C  to 
500  C  was  chosen.  All  furnace  temperatures  mentioned  were 
maintained  within  five  degrees  of  the  temperature  stated. 

The  heat  treatments  for  the  T  series  were: 

1.  A  number  of  samples  were  sealed  together  in  an  evacuated 
quartz  ampoule,  solution  heat  treated  at  1050  C  for 

24  hours,  and  cold  water  quenched,  breaking  the  quartz 
ampoule  under  water. 

2.  The  samples  were  then  individually  sealed  in  evacuated 
quartz  ampoules  aged  as  indicated  in  Table  3.2,  and 
quenched  in  cold  water. 

The  length  of  ageing  was  determined  by  assuming  that  the 
square  of  distance  the  solute  travels  is  approximately  equal 
to  the  diffusion  coefficient  for  the  solute  times  the  time 
available.  In  all  cases  a  phosphorus  atom  was  allowed  on 
average  to  move  at  least  one  micron,  based  on  a  low 
temperature  extrapolation  of  a  high  temperature  diffusion 
coefficient  [32].  The  paramagnetic  diffusion  coefficient  may 
be  used  for  higher  temperatures,  but  extremely  long  times 
are  necessary  at  low  temperatures  due  to  a  diminished 
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Table 


.2  Ageing 

temperature  and  times 

for  final  samples 

Sample 

Ageing  Time 

Diffusion  Distance 

( hrs ) 

(microns ) 

T800C 

240 

43.0 

T750C 

504 

33.4 

T700C 

960 

23.1 

T650C 

1896 

15.1 

T600C 

3840 

9.1 

T550C 

4272 

3.7 

T500C 

4920 

1.3 
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the  Curie  temperature  (see  section  4.11). 


3.3  Metal lographic  Specimen  Preparation,  Examination,  and 
Chemical  Analysis 

All  samples  (buttons)  were  sectioned  using  a  Buehler 
Low  Speed  Saw  with  a  diamond  blade.  Examples  of  buttons  are 
shown  in  Plate  3.6.  The  sectioned  samples  were  then  mounted 
in  a  1  in.  mold  using  a  copper - f i 1 1 ed  thermo-setting  epoxy; 
this  allowed  the  mount  to  conduct  electrons,  a  useful 
property  when  viewing  with  a  Scanning  Electron  Microscope 
(SEM) . 

After  mounting,  all  samples  were  polished  in  the 
following  sequence  using  a  Buehler  Minimet  Automatic 
Pol i sher : 

1.  420  grit  silicon  carbide  paper  and  water, 

2.  600  grit  silicon  carbide  paper  and  water, 

3.  6  micron  diamond  paste  and  oil  on  Buehler  AB  Texmet 

paper , 

4.  3  micron  diamond  paste  and  oil  on  Buehler  AB  Texmet 
paper , 

5.  1  micron  diamond  paste  and  oil  on  Buehler  AB  Texmet 
paper . 

A  final  polish  with  Linde  B  alumina  was  attempted  but  the 
specimens  were  found  to  pit  excessively.  All  specimens  were 
given  a  final  10  second  etch  with  a  2%  Nital  solution. 
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Plate  3.6  An  example  of  final  alloy  (buttons)  produced 
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Micrographs  were  taken  of  the  sample  mi crostructures  at 
various  magnifications  using  a  Zeiss  Ultraphot  III 
Meta  1 lograph .  The  samples  were  also  viewed  with  much  greater 
resolution  using  an  ISI  60  Scanning  Electron  Microscope 
(SEM).  This  SEM  is  fitted  with  a  National  Semiconductor 
Energy  Dispersive  Spectrometer  (EDS)  and  a  Tracor  Northern 
T N - 17 10  multichannel  analyser.  These  instruments  produce  an 
X-ray  spectrum  of  the  elements  present  under  the  electron 
beam,  permitting  a  qualitative  chemical  analysis. 

Following  are  a  few  examples  of  the  specimen 
microstructures  and  X-ray  spectra.  Plate  3.7  shows  the 
microstructure  of  the  produced  iron  phosphide  after  vacuum 
arc-melting  into  a  button.  The  nominal  20  wt%  P  present 
would  place  this  alloy  in  the  Fe3P-Fe2P  range.  Grains  with 
higher  P  content  are  char acter i zed  by  their  rectangular 
morphology.  Plate  3.8  shows  a  typical  X-ray  spectrum  from 
this  phosphide.  As  is  evident,  no  impurities  are  detectable. 
Plate  3.9  shows  a  representative  microstructure  of  T700C. 

The  phosphide  precipitates  appear  as  plates  and  needles 
along  specific  crystallographic  planes  and  directions  [24]. 
Also  present  are  grain  boundary  a  1 lotr i omorphs  as  described 
by  Aaronson  and  co-workers  [33,34].  Plate  3.10  shows  a 
representative  X-ray  spectrum  from  the  bulk  solid  solution 
of  T700C.  The  large  peak  is  an  iron^-alpha  peak.  Other  than 
iron  and  phosphorus  there  are  no  detectable  elements. 

Plate  3.11  shows  a  representative  microstructure  of 
T600C .  This  plate  is  similar  to  Plate  3.9,  with  all  relevant 
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Plate  3.7  Microstructure  of  produced  iron  phosphide,  X520 
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Plate  3.8  A  typical  X-ray  spectrum  for  the  produced  iron 
phosphide 


Plate  3.9  Precipitation  in  T700C, 


X520 
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Plate  3.10  A  typical  X-ray  spectrum  from  the  bulk  solid 
solution  of  T700C 


> 
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Plate  3.11  Precipitation  in  T600C,  X480 
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features  the  same.  However,  the  precipitates  are  smaller  and 
less  widely  spaced.  This  is  as  expected  since  the  diffusion 
coefficient  of  P  decreased  more  rapidly  than  the  time 
available  for  diffusion  increased.  Plate  3.12  shows  the 
mi crostructure  of  sample  T500C.  Here  the  precipitates  have 
had  insufficient  time  to  grow  therefore  they  are  very  small 
and  finely  spaced.  This  may  be  explained  by  considering  the 
very  recent  work  of  Luckman  et  al.  [35]  which  shows  that  P 
diffusivity  decreases  more  rapidly  than  predicted  from  high 
temperature  extrapolation. 

A  quantitative  chemical  analysis  of  the  specimens  was 
carried  out  on  an  ARL  electron  microprobe  in  the  Department 
of  Geology,  by  courtesy  of  Dr.  D.  G.  W.  Smith.  Established 
experimental  techniques,  technical  assistance  and  data 
analyses  were  provided.  An  EDS  system  was  first  used  to 
determine  the  purity  of  the  samples.  Then,  because  only  the 
two  elements  needed  analysis,  a  Wavelength  Dispersive 
Spectrometer  (WDS)  was  used  for  all  quantitative  analyses. 
The  resulting  raw  data  (X-ray  counts  for  1C -alpha  peak  and 
background  for  each  element)  were  analysed  by  an  iterative 
Atomic  Number-Atomic  Absorpt i on- F 1 uorescence  correction 
program  called  EDATA2.  Sensitivity  of  the  system  was  better 
than  0.02  at%  of  P. 

The  experimental  results  are  given  in  section  4.4.2 
where  they  are  used.  By  comparing  100  random  areas,  the 
produced  master  alloy  phosphides  were  analysed  as  64%  Fe2P 
and  36%  Fe^P  to  give  an  average  18.6  wt%  P  (a  loss  of 
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Plate  3.12  Precipitation  in  T500C, 


X480 
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1.4  wt%  P  from  the  nominal  20  wt%  P  in  production).  This 
value  was  used  to  produce  the  nominal  1.5  wt%  P  final  alloys 
(~1.30  wt%  P  actually  present).  All  final  alloys  were  vacuum 
arc-melted  twice  to  assure  homogeneous  distribution  of  P. 
This,  plus  the  fact  that  any  accidental  overheating  of  the 
melt  increases  its  vapour  pressure,  accounts  for  the  loss  of 
phosphorus . 


3.4  Theoretical  Techniques 

The  major  part  of  this  thesis,  presented  in  chapter  4, 
is  a  theoretical  analysis  of  i ntergranul ar  segregation  in 
iron.  The  equipment  and  techniques  used  are  described  here. 

3.4.1  Equipment 

All  calculations  were  executed  on  a  Hewlett  Packard 
(HP)  9825A  Computer  using  the  HP  Basic  programing  language. 
Most  graphs  were  produced  on  an  HP  9872A  Graphic  Plotter. 

The  programs  used  were  monitored  using  a  HP  9871A  Line 
Pr i nter  . 

3.4.2  Numerical  Analysis 

All  analytical  equations  were  solved  directly. 
Generally,  these  functions  (chemical  potentials,  limit  of 
solubilities  and  Case  I  segregation  equations  of  chapter  4) 
were  plotted  on  graphs  immediately  as  they  were  calculated. 
However,  there  was  also  the  need  to  simultaneously  solve  the 
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roots  of  two  or  more  equations  (e.g.  Case  II  magnetic 
segregation  free  energies,  segregation  curves  and  boundary 
Curie  temperatures  of  chapter  4).  The  finding  of  roots  of 
equations  of  the  form  f(x)  =  0  is  one  of  the  most  frequently 
occurring  problems  in  scientific  work.  In  general,  it  is 
only  possible  to  obtain  approximate  solutions,  relying  on 
some  computation  technique  to  produce  the  approximation. 

A  linear  iteration  method  was  employed  for  finding  the 
roots  of  simultaneous  equations.  For  a  single  equation, 
f ( x )  =  0  may  be  expressed  as  x  =  g(x)  and,  geometrically,  a 

root  of  f(x)  =  0  is  a  number  x  =  z  for  which  the  line  y  =  x 

intersects  the  curve  y  =  g(x).  The  equation  x  =  g(x)  is 
solved  by  a  converging  recursion  of  x.  It  must  be  noted  that 
the  condition  |g'  (x)|  <  1  should  be  observed  for 
convergence . 

For  simultaneous  equations,  consider  the  following 
example  with  2  equations: 

f ( x , y )  =  0  and  x  =  F ( x , y ) 

g(x,y)=  =  0  and  y  =  G(x,y). 

Sufficient  but  not  necessary  conditions  [36]  for  the 
convergence  of  the  linear  iteration  method  are: 

1.  F  and  G  and  their  first  derivatives  are  continuous  in  a 
neighbourhood  R  of  the  desired  root. 

2.  The  sum  of  the  absolute  values  of  the  partial 
derivatives  of  F  and  G  must  be  less  than  one  for  all 
points  (x,y)  in  R. 

3.  The  initial  approximation  of  (x,y)  is  chosen  in  R. 
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Whenever  possible  an  accelerated  technique  was  used  in  place 
of  the  linear  iteration  method.  This  technique  is  called  the 
Wegstein  method  for  algebraic  convergence  [37]  and  is  based 
on  a  projection  technique.  That  is,  a  projection  from  two 
Known  points  on  the  x  =  f(x)  line  determines  a  point  on  the 
y  =  x  line.  This  x  value  is  then  used  for  f(x).  The 
projection  requires  the  solution  of 

x  -  -  y]x2 _ 

3  X1  "  x2  +  y2  "  yl 

and  x^  and  x^  now  become  the  new  Known  points.  This 
procedure  is  repeated  until  a  rapid  convergence  is  obtained. 
The  numerical  problems  encountered  in  section  4  required  the 
solving  of  at  least  two  simultaneous  equations  such  as: 

x  =  f  (  x  ,  y ,  T  ) 
y  =  g  (  x  ,  y ,  T  ) 

This  was  accomplished  by  up  to  four  independent  iteration 
techniques  which  will  be  called  A,  B,  C,  and  D: 

A.  For  each  constant  T  required,  a  linear  iteration 

technique  for  finding  x  and  y  was  used.  If  A  failed  (a 
certain  number  of  recursions  was  exceeded),  control  was 
passed  to  B.  This  usually  happened  when  the  slope  of  the 
curve  became  very  steep. 

B.  One  of  the  curves  was  'inverted'  and  the  following 

equation  solved: 

T  =  h  (  x  ,  y ,  T  ) 
y  =  g(x,y,T) 


46 


For  each  constant  x,  a  linear  iteration  technique  for 
finding  y  and  T  was  used.  If  B  failed,  control  was 
passed  to  A.  However,  if  A  and  B  both  failed,  control 
was  passed  to  C.  A  and  B  both  failing  occurred  often  and 
usually  indicated  an  oscillation  in  the  iteration. 

C.  This  iteration  was  useful  only  when  a  special  condition 

applied,  e.g.  when  x  +  y  =  1.  For  a  constant  x,  y  was 
set  and  T  could  be  found  using  a  Wegstein  acceleration. 
If  C  could  not  be  used,  control  was  passed  to  D. 

D.  This  technique  consisted  of  a  Wegstein  accelerator  inner 

iteration  and  a  graphical  bisection  outer  iteration. 

That  is,  for  each  step  of  x,  an  iteration  similar  to  B 
found  T  and  y.  The  value  of  x  was  stepped  until  a  change 
in  sign  for  the  curves  was  noticed.  Each  step  of  x  was 
then  halved  as  the  value  of  x  oscillated  around  the  root 
until  a  certain  desired  accuracy  was  reached.  This 
technique  was  very  slow  but  always  produced  a  root.  If  a 
special  condition  occurred,  then  control  would  be  passed 
to  C. 

The  segregation  and  free  energy  curves  of  section  4  were 
produced  by  Program  1,  Appendix  2.  A  flow  chart  of  the 
program  is  also  presented  in  Appendix  2.  When  an  iterative 
method  was  used,  up  to  240  points  were  calculated  per 
function  to  produce  smooth  curves.  This  data  was  stored  on 
magnetic  tape  during  calculation  and  later  retrieved  for 
plotting. 


4.  RESULTS  AND  DISCUSSION 
This  chapter  contains  a  theoretical  thermodynamic 
development,  some  sections  of  which  have  been  previously 
published  by  the  author  [38].  The  effects  of  ferromagnetism 
on  selected  properties  of  iron-based  solid  solutions  are 
approached  from  a  macroscopic  equilibrium  thermodynamic 
viewpoint.  These  properties  include  the  solid  solubility 
limit  of  a  solute,  the  intergranular  segregation  of  a 
solute,  and  diffusion  within  the  solution.  Although  the 
model  described  is  applicable  to  any  solute  in  iron, 
phosphorus  in  iron  is  used  as  an  example  and  the 
experimental  results  for  the  behavior  of  phosphorus  in  alpha 
iron  are  presented  where  convenient  to  support  the 
thermodynami c  development.  Starting  from  basic  definitions, 
the  model  will  be  built  up  to  deal  with  binary  and  ternary 
solid  solutions.  Although  all  symbols  used  are  defined  at 
the  time  of  their  first  appearance  in  the  text,  a  complete 
compilation  of  symbols  and  their  definitions  is  given  in 
Appendix  1.  Although  use  was  made  of  the  work  of  others, 
particularly  Nishizawa  and  Guttmann,  the  treatment  is 
original  unless  otherwise  stated. 


4.1  The  Free  Energy  Function  for  alpha  Iron 

The  metal  under  consideration  is  body  centered  cubic 
alpha  iron.  In  the  pure  state  alpha  iron  is  stable  from 
0  to  1185  K.  When  dealing  with  conditions  of  equilibrium,  it 
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is  convenient  to  describe  the  system  in  terms  of  state 
variables  or  functions.  As  the  state  of  the  system  is 
changed,  the  change  of  any  state  function  depends  only  on 
the  value  of  the  function  in  the  initial  and  final  states, 
i .e.  the  change  in  the  functions  is  independent  of  the  path 
along  which  the  change  has  occurred.  State  variables  are  not 
necessarily  independent  of  one  another.  However,  once  the 
state  of  a  system  is  specified  by  the  values  of  a  few  state 
functions,  the  values  of  all  are  fixed. 

Consider  the  Gibbs  free  energy  function  for  a  pure 
substance.  The  Gibbs  free  energy,  G,  is  a  state  function 
defined  by 


Entropy,  S,  is  a  state  function  and  enthalpy,  H,  is  defined 
by 

H  =  U  +  PV 

(4.2) 

where  internal  energy,  U,  pressure,  P,  and  volume,  V  are 
state  functions.  For  a  pure  substance,  G  depends  on  the 
chemical  nature  of  the  substance,  the  amount  of  the 
substance  in  the  sample,  the  state  (solid,  liquid,  gas)  in 
which  it  exists,  the  temperature  and  applied  pressure.  The 
free  energy  of  all  substances  decreases  with  increase  in 
temperature.  At  0  K  the  free  energy  is  equal  to  the  total 
internal  energy  of  a  solid,  which  is  equivalent  to  its 
enthalpy.  An  increase  in  temperature  increases  the  enthalpy. 
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The  rate  of  increase  of  the  enthalpy  with  temperature  is  the 
heat  capacity  of  the  substance.  Similarily,  as  free  energy 
decreases  with  increasing  temperature  the  rate  of  decrease 
is  the  entropy  of  the  substance.  In  general,  the  Gibbs  free 
energy  may  be  given  by 


G  =  °G  +  RT£na 


(4.3) 


where  °G  is  the  free  energy  of  the  substance  when  present  in 
its  standard  state  and  'a'  is  its  activity. 

Consider  the  total  Gibbs  free  energy  function,  G^6,  for 
alpha  iron.  By  definition, 


GFe  =  HFe  -  TSFe 


(4.4) 


However,  as  pointed  out  by  Zener  [39]  and  Weiss  and  Tauer 
[40],  the  total  free  energy  function  should  be  analysed  by 
its  component  parts,  i.e. 


GFe  =  G Fe  .  +  GFe...  +  GF<f 

magnetic  lattice  electronic  ^ 

This  type  of  approach  was  first  taken  by  Zener  [39]  who 
explained  the  abnormal  form  of  the  gamma- loop  in  the  Fe-Cr 
system  as  being  due  to  the  magnetic  transition  in  alpha 
iron.  Later,  Weiss  and  Tauer  [40,41]  numerically  evaluated 
the  magnetic  (subscript  mag)  and  non-magnetic  (subscript  nm) 
components  of  enthalpy  and  entropy  (and  hence  of  free 
energy)  for  alpha  iron.  Both  the  enthalpy  and  the  entropy 
may  be  estimated  from  the  heat  capacity  (or  specific  heat). 
Although  below  the  Curie  temperature  the  electronic  heat 
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capacity  for  each  direction  of  electron  spin  is  different 
and  should  be  considered  separately,  the  overall  lattice  and 
electronic  heat  capacities  are  similar.  Therefore,  the  total 
Gibbs  free  energy  may  be  given  by 


,Fe 

u 


G 


Fe 

nm 


+ 


G 


Fe 
ma  a 


and  because 


.Fe 


1  s 


an  extensive  function, 


(4.6) 


QFe  =  H Fe  _  T$ Fe 

mag  mag  mag  (47) 

The  term  magnetic  enthalpy,  as  an  example,  refers  to  the 
magnetic  component  of  enthalpy.  There  is  an  excellent  review 
by  Miodownik  [42]  which  deals  with  non-linear  magnetic 
components  of  free  energy  and  phase  stability  in  magnetic 
a  1 loys . 

The  magnetic  enthalpy  and  entropy  functions  used  in  the 
present  work  are  based  on  the  data  of  Weiss  and  Tauer  [40] 
as  shown  in  Table  4.1. 

From  this  data,  empirical  functions  of  magnetic  enthalpy  and 
entropy  were  developed  by  the  author  as  described  in  the 
following  paragraphs. 

Phase  transitions  are  generally  classified  by  order 
according  to  the  presence  of  a  discontinuity  in  that  order 
of  derivative  of  the  free  energy  with  respect  to  a  given 
state  variable  such  as  temperature.  First  order  transitions 
(a  discontinuity  is  present  in  the  first  derivative)  have 
discontinuities  in  entropy,  volume  and  enthalpy  at  the 
transition  temperature.  Examples  of  first-order  transitions 
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Table  4.1  The  magnetic  thermodynamic  functions  of  alpha  iron 
[40] 


Temperature 

Enthalpy 

Entropy 

Free  Energy 

(K) 

( J/mole) 

( J/mole/K ) 

( J/mole ) 

0 

-8732 

0.00 

-8732 

400 

-8636 

0.21 

-8732 

500 

-8535 

0.38 

-8720 

600 

-8368 

0.67 

-8725 

700 

-8037 

1  .  17 

-8770 

800 

-7518 

1.88 

-8856 

900 

-6538 

3.01 

-9022 

1000 

-4831 

4.86 

-9247 

1100 

-2671 

6.99 

-9691 

1200 

-1842 

7.70 

-10360 

1300 

-1398 

8.  12 

-11082 

1400 

-1122 

8.33 

-11954 

1500 

-925 

8.46 

-12784 

1600 

-787 

8.54 

-13615 

1700 

-716 

8.58 

-14451 

1800 

-674 

8.62 

-16190 

52 


are  fusion,  vaporization  and  allotropic  transformations. 
Second-order  transitions  have  discontinuities  in  the  heat 
capacity,  the  compressibility  and  the  thermal  expansion 
coefficient.  Some  examples  of  second-order  transformations 
are  order -di sorder  transformations  in  alloys  and  the  onset 
of  ferromagnetism. 

For  materials  with  a  paramagnetic  to  ferromagnetic 
transformat  ion  (e.g.  alpha  iron)  the  magnetic  heat  capacity 
has  a  very  large  peak  (spike)  at  the  Curie  temperature  "Te 
as  shown  in  Figure  4.1.  Some  previous  work  has  simulated  the 
magnetic  heat  capacity  Cp^  [43]  using  fitting  formulae 
such  as 


p  —  _  /  r  t  i  <  n,i  ,  r  i  i  n  t  \ 

CPmag  -  Arv£rt{  LI  +  L  J/Ll-L  j) 


=  BRc.:([l+L‘m]/[l-L‘"']),  L  >  1 


-m 


(4.8) 


where  L  =  T /fe  and  n  =  3,  m  =  5.  However,  in  the  present 

c 

work  these  formulae  were  fitted  to  magnetic  heat  capacity, 
enthalpy  and  entropy  data  through  the  following: 


Fe 

fcoFe 

dT 

mao 

J  mag 

Fe  = 

fCDFe 

dT/7 

mac 

J  mag 

Since  the  heat  capacity  of  iron  is  infinite  at  the  Curie 
temperature  an  analytic  determination  of  the  integrals  in 
equations  (4.9)  is  not  possible.  Instead,  a  power  series 
expansion  for  the  natural  logarithmic  terms  in  equations 


(4.8)  was  used,  where 
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CM/3H0 W/r>  AlIGVdVD  ±V3H  DI13NDVW 


F igune  4 . 1 

The  magnetic  heat  capacity  of  alpha  iron. 
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£h( [1 +X ]/ [ 1 -X] )  =  2[X+X3+X5+ . ] 

I  5 

(4.10) 

This  allowed  a  simple  integration  of  equations  (4.9)  and  a 
fitting  of  equations  (4.9)  to  the  data  in  Table  4.1  with  the 
following  results:  for  L  <  1,  A  =  1.08444  and  n  =  3.613  and 
for  L  >  1,  B  =  0.67423  and  m  =  4.766.  This  calculated 
magnetic  heat  capacity  of  alpha  iron  is  shown  in  Figure  4.1. 
The  discontinuity  at  the  Curie  Temperature  was  avoided  by 
arbitrarily  setting  a  maximum  for  the  heat  capacity. 

Figure  4.2  shows  the  magnetic  entropy  curve  while  Figure  4.3 
shows  the  magnetic  enthalpy  and  resultant  Gibbs  free  energy 
curves  for  alpha  iron  (data  from  Table  4.1).  An  assumed 
maximum  value  for  magnetic  heat  capacity  at  the  Curie 
temperature  corresponds  to  straight  line  approximations 
immediately  about  the  Curie  temperature. 

The  magnetic  free  energy  of  alpha  iron,  C-j^®  (T)f  as 
shown  in  Figure  4.3,  may  be  divided  by  the  Curie  temperature 
into  two  regions,  an  uncoupled  (paramagnetic)  state,  u  gF®  , 
and  a  coupled  ( ferromagnet  ic )  state,  G>  ,  such  that,  in 


mag 


approx i mat  ion : 


iFe 

'mag 


-AT, 

T  > 

TFe 

c 

-ATFe , 

T  < 

TPe 

(4.11) 
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Figure  4.2 

The  magnetic  entropy  of  alpha  iron  (+  denotes  data  [40]) 
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Figure  4.3 

The  magnetic  enthalpy  and  Gibbs  free  energy  of  alpha  iron  (+  denotes  data  [40]). 
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The  constant  A6  may  be  evaluated  by  considering  equation 
(4.7),  and,  differentiating  with  respect  to  T  when  L  >  1 
gives 


9G 

IT 


Fe 

mag 


9H  Fe  -  T9S  Fe  -  SFe 
mag  mag  mag 

9  T  9T 


A 


Evaluating  the  components  at  L  >  1  shows 

i  Fe 


mag 

Fe 

mag 


-*  0  J/mole 

=  9  J/mole/K 


U$Fe 

mag 


which  implies  that 

SF^g(T>>TFe)  ->  A  ,  a  constant 

uGFe  =  -TuSFe 
mag  mag 

crFe  _  T Fe  ucFe 
mag  c  mag 

An  alpha  Fe-X  solid  solution  will  now  be  considered. 


(4.12) 


(4.13) 


(4.14) 


4.2  The  alpha  Fe-X  Solid  Solution 

By  analogy,  the  free  energy  of  an  alpha  Fe-X  solid 
solution7,  Ga ,  is  assumed  to  have  magnetic  (subscript  mag) 
and  non-magnetic  (subscript  nm)  components  [39,40,44],  such 
that 


6This  is  not  the  same  constant  as  in  equation  (4.9). 

7The  solution  has  mole  fractions  X  and  X  of  the  components 
i  and  j . 


' 
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Ga  =  Ga  +  Ga 

nm  mag  ^ 

Consider  the  formation  of  a  binary  solution  from  the  pure 
components  in  their  standard  state.  The  reaction  is 


X.i  +  X  j  =  (X  ,X  )ij  +  AG 
l  J  i  J  m  (4.16) 

With  this  reaction  is  associated  a  free  energy  change  called 

the  free  energy  of  mixing,  AGm ,  such  that 


AG  =  Ga  -  X  .'V  -  X  .°GJ 
m  l  i 


(4.17) 


For  a  solid  solution 


Ga  =  y 1  X  .  +  X  . 

1  J  (4.18) 

where  jj1  is  the  chemical  potential  of  component  i. 

Therefore,  the  free  energy  of  mixing  is 

AGm  =  X .  [u1  -  V]  +  Xj [uJ  -  °Gj]  (4.19) 


It  is  assumed  that  the  non-magnet ic  component  of  the 
free  energy  of  the  solid  solution  is  approximated  by  the 
regular  solution  model.  Regular  solutions  are  solutions 
formed  for  components  for  which  the  entropy  of  mixing  has 
the  value  expected  for  an  ideal  solution  but  for  which  the 
enthalpy  of  mixing  is  not  zero.8  A  non-zero  enthalpy  may  be 
expressed  by  a  power  series  of  the  component  concentrations. 


8The  enthalpy  of  mixing  is  zero  for  an  ideal  solution. 
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Therefore,  as  expressed  by  Nishizawa  et  al.,  the 
non-magnet ic  component  of  free  energy  is 


g“  =  [l-Xy]V®  +  X  °GX  +  Xv[l-Xv]IFeX 
nm  X  nm  X  nm  XL  XJ  nm 

+  RT([1-X  ]£w(l-Xx)  +  Xx£nXx> 


(4.20) 


-This  free  energy  is  expressed  as  an  atomic  fraction  of 
solute  X,  Xx,  where  °G^m  is  the  standard  free  energy  for 
non-magnetic  i  and  ijjj  is  the  non-magnetic  interaction 
parameter  between  components  i  and  j. 

Unfortunately,  it  is  not  possible  to  theoretically 
estimate  the  magnetic  free  energy  in  the  same  way;  instead, 
descriptive  modifications  to  the  magnetic  free  energy  of 
alpha. iron  are  made  [40,45],  as  described  now. 

It  is  well  Known  that  the  addition  of  many  species  of 
solute  to  alpha  iron  alters  its  Curie  temperature.  This  may 
be  observed  in  any  major  compilation  of  binary  phase 
diagrams  such  as  found  in  the  Metals  Handbook  [47].  The 
change  in  Curie  temperature  of  a  solution  with  change  in 
solute  content  is  often  small  and  consequently  easy  to 
ignore.  To  ignore  a  shift  in  the  Curie  temperature  is,  in 
effect,  making  a  zeroth  approximation. 

Nishizawa  et  al.  [44]  considered  that  the  Curie 
temperature  of  iron  is  shifted  from  1043  K  by  dissolving  the 
atomic  fraction  X  of  solute  X,  such  that 


Ta  =  jFe  +  ATYX 
C  C  XX 


(4.21 ) 


The  magnetic  parameter  for  solute  X,  aTx  ,  represents  the 
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effectiveness  of  a  given  solute  in  changing  the  Curie 
temperature.  This  linear  approximation  would  be  valid  at 
least  for  dilute  solid  solutions.  Assuming  such  a 
concentration  dependence  for  the  Curie  temperature  Hillert 
et  al.  [45,46]  have  shown  that  the  magnetic  free  energy  for 
alpha  Fe-X  may  be  obtained  by  modifying  the  free  energy  of 
alpha  iron.  Figure  4.4  shows  the  modification  where  the 
magnetic  free  energy  curve  for  alpha  Fe-X  is  displaced  along 
the  paramagnetic  free  energy  line,  U(^g  •  That  is,  the  alpha 
Fe-X  magnetic  free  energy  curve  maintains  the  same  position 
with  respect  to  the  new  Curie  temperature  as  the  alpha  iron 
free  energy  curve  to  the  old  Curie  temperature  [44]. 

The  displacements  mentioned  may  be  described  by 
considering  for  alpha  iron  a  temperature  T-j  such  that 

GFe  (T, )  =  -TFe  uSFe 
magi  c  mag 

Fe 

where  T  -  T,  =  K,  a  constant.  . . 

c  1  (4.22) 

For  the  solid  solution  alpha  Fe-X,  there  exists  a 

temperature  T^  such  that 


where 

This  implies  that 


■  a 

'mag 


(TJ 


Ta  tuFe 
c  maq 


K 


Ta  =  T^e  _  T  +  T 

c  c  1  l2 
Ti  *  T2  -  ATXXX 


(4.23) 


(4.24) 


For  alpha  Fe-X,  T£  is  the  system  temperature,  T,  but  if 
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Ta, 


and 


then  (see  Figure  4.4) 

T*  =  T  -  ATXXX 

Ga  (T)  =  G Fe  ( Ta ) 
magv  '  magu  ' 


ATXXX 


uSFe 

mag 


(4.25) 


The  above  is  true  for  all  values  of  T. 

Weiss  and  Tauer  [41]  pointed  out  the  influence  of  the 
solute  upon  the  size  of  the  magnetic  component,  introducing 
what  may  be  called  a  solute  intensity  factor  M  ,  so  that  in 

A 

an  alpha  Fe-X  solid  solution 


Graag  -  n-MYX¥][G^„(T°)  -  ATvXvuS„Ff J 


rx 


mag 


X  X  mag J  (4.26) 


Therefore,  the  free  energy  of  the  alpha  Fe-X  solid  solution 
may  be  expressed  as 

Ga  =  [l-Xv]°GFf  +  X„°GX_  +  Xv[l-Xv]IFeX 


nm 


nm 


nm 


+  RT{ [ 1 -XY ]£n( 1 -X v )  +  Xv  £nXv} 


+  [1-MXX  ][GF®  (Ta)  -  ATyX  US Fe  ] 
a  a  mag  X  X  magJ 


(4.27) 


In  recent  work  by  Nishizawa  et  al.  [44],  the  interaction 
term  has  been  reformulated  to  include  a  magnetic 
interaction.  In  regard  to  a  dilute  solid,  it  was  pointed  out 


that : 
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IFex 

TFeX 
x  nm 

+  iFeX 
mag 

jFeX 

mag 

ATX[S 

Fe  (Ta}  .  u  Fe 
mag  '  magJ 

c .  FeX  — 
mag 

j  FeX 

1  nm 

-  ATyuSFe 

X  mag 

u  j  FeX 

1  mag 

j  FeX 

1  nm 

(4.28) 


The  above  was  justified  in  the  following  manner  by 
considering  Hillert  et  al.  [45,46].  The  excess  free  energies 
or  the  interaction  free  energies,  I,  are  very  small  at  low 

alloy  contents.  For  the  magnetic  term,  GFe  ,  a  small  alloy 

mag 

addition  AX  causes  the  Curie  temperature,  Tc  ,  to  be 
displaced  by  aT  .  Therefore,  the  change  in  free  energy  of 
the  solution  at  any  temperature  (per  mole  alpha  iron)  is 


AT  3GFe 
c _ mag 

3T 


which,  per  mole  alloying  element,  is 


lTc  <ag 
AX  3  T 


In  the  present  case 


AX  =  X 


AT  =  - AT v  X 


XAX 


and  term  (4.30)  becomes 


ATXS 


Fe 

mag 


(T) 


(4.29) 


(4.30) 


(4.31 ) 


(4.32) 


Now,  considering  equation  (4.26), 
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jFeX 

mag 


ATY[SFe 
X  mag 


uSFe  ] 
mag J 


(4.33) 


which  is  the  second  of  Nishizawa's  equations  (4.28). 
Therefore,  the  approximated  alpha  Fe-X  free  energy  [44]  may 
written  as 


i°  =  [1  -Xv  ]°GFe  +  Xy°GX  +  IFeXXv[l-Xv] 
L  X  nm  X  nm  mti  X  X 


+  RT{ [ 1  - X^ ]£«( 1  - X  )  +  Xx£nXx} 

+  ^VxKag'1’  +  ^Vx^g^ 


u  r  Fe  -|  -i 

mag  J  (4.34) 


where  GFe  (Ta)  =  GFe  (T)  +  S F®  ( T01) £TyX v 

mag  mag  mag'  X  X  (4.35 

The  above  free  energy  may  be  used  to  determine  the  chemical 
potentials  of  the  constituents  of  the  solid  solution. 


4.3  Chemical  Potentials  of  alpha  Fe-X  Solid  Solutions 

The  free  energy  of  an  alpha  Fe-X  solid  solution  may  be 
expressed  in  terms  of  its  component  chemical  potentials  as 
expressed  in  equation  (4.18).  Taking  the  differential  of 
equation  (4.18)  gives 

dGa  =  yFedXpe  +  yXdXx 

where  Xc  dyFe  +  XvdyX  =  0  (Gibbs-Duhem  equation), 

• e  a  (4.36) 

In  a  binary  solution  dXpe  =  -dXx  and  therefore  the  chemical 

potentials  are  described  by 
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UtQ  =  Ga  -  EX.^GJ'  =  Ga  -  XY  3Ga 
j  J^Xj  A  3  X  x 

yX  =  Ga  -  ZX  3Gj  +  iG01  =  uFe  aGa 

j  SXj  3XX  3XX  (4.37) 

These  chemical  potentials  may  be  separated  into  magnetic  and 
non-magnetic  components.  The  partial  differential  of  the 
free  energy  of  the  solution  with  respect  to  the  solute 
content  is 


3  G 
3  X 


-°GPe  +  [i-Xy]3°GFe  +  °GX  +  X  Y  3 0  G X 
nm  L  XJ  r.m  nm  X  nrn 


3  X 


3  X 


+  [l-2Xv]I^  +  Xv[i-x  ]3T=  +  RT[drtX  y-£n(l-X  )] 


X J  nm 


X 


X J.  nm 
3  X 


X 


'  £VXX3MX][Gmag<Ta>  ' 


3 


+  n-MxXx][<^(T“)  -  <S^g] 


3  X 


X 


(4.38) 


At  this  point,  it  is  necessary  to  evaluate  the  terms 
present : 

1.  The  standard  non-magnetic  free  energy  of  pure  alpha  iron 
is  a  function  of  temperature  only,  such  that, 

3°GF*  ■  0 
nm 

3  X  x  (4.39) 

2.  Also,  the  standard  non-magnetic  free  energy  of  pure  X  is 
a  function  of  temperature  only,  such  that, 
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a°Gx  =  o 
nm 

3XX  (4.40) 

3.  The  non-magnetic  interaction  parameter  is  a  function  of 
temperature  only,  such  that 

aiFeX  =  o 

nm 

3XX  (4.41) 

4.  The  modified  magnetic  free  energy  of  alpha  iron  is  a 
function  of  temperature  and  the  component 
concentrations,  such  that 

3GFe  ( Ta)  =  3GFe  ( Ta )  3Ta 
mag v  '  _ mag  v  ‘  %y 

3  X„  ,Ta  d*X 


3T 


where 


•  a 


3T  =  -AT 
3  X  v 


(4.42) 


The  modified  magnetic  free  energy  is  defined  as  a 
function  of  modified  temperature  by  reference  to 
equation  (4.7).  Differentiating  equation  (4.7)  with 
respect  to  the  modified  temperature  gives 

3GFe  (Ta)  =  3HF®  (Ta)  -  SF®  (Ta)  -  Ta3SFe  (Ta) 
mag v  '  _ mag  magv  mag  ' 

3Ta  3Ta  3Ta  (4.43) 

As  stated  previously,  the  magnetic  enthalpy  and  entropy 

were  obtained  from  the  magnetic  heat  capacity  in  such  a 


manner  that 
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3HFe  (Ta)  =  CpFe  (Ta) 
mag  v  '  magv  ' 


9T 


a 


3SF®  (Ta)  =  CpF®  (Ta) 
mag  mag 


3T 


a 


■  ct 


3GFe  ( Ta)  =  ATvSFe  (Ta) 
mag v  '  X  magv 

3XX  (4.44) 

5.  The  modified  magnetic  entropy  may  be  handled  in  a 
similar  way,  such  that 


!fiaiTa)  =  -ATxCp^q(Ta) 

3XX  Ta  (4.45) 

However,  as  noticed  in  figure  4.1,  the  modified  heat 
capacity  for  alpha  iron  has  small  values  even  near  the 
transition  temperature.  Therefore,  this  term  has  been 
ignored  throughout. 

6.  The  solute  intensity  factor  may  be  a  function  of 

temperature  and  solute  concentrations.  In  this  case, 
was  chosen  to  be  a  constant,  M,  which  depends  on  the 
magnetic  properties  of  the  solute  as  discussed  below. 

Equation  (4.38)  may  now  be  written  as 


3  Ga 
3  X  X 


-°GFe  + 
nm 


or  X 
Gnm 


+  [1-2Xx]IF®X+  RTU>.XX-M1-XX)] 


-  -  ATXXX  USnP 

+  l>'VAVSIag<T“>  -  “#g] 


(4.46) 


The  validity  of  the  approximation  in  equation  (4.36)  may  be 
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seen  from  the  fact  that  any  derivative  of  the  solution  free 
energy  remains  unchanged  with  the  approximation. 

The  chemical  potentials  of  both  Fe  and  X  in  the  solid 
solution  can  also  be  separated  into  magnetic  and 
non-magnetic  components.  Therefore,  the  chemical  potential 
of  Fe  in  alpha  Fe-X  is 


F  e 

Fe  ,  Fe 

u  +  u 

nm  mag 

Fe  = 

0GFe  +  Xyl 

nm 

nm  X 

Fe 

=  GFe  ( Ta ) 
mag  ' 

mag 

FeX 

nm 


+  RUn(l- 


at  y  ucFe 

“  ATyXy  O 

X  X  m  a  9 


-  Xx[l-MXx]ATx[S^g(T“) 


usFe  ] 

mag  J 


(4.47) 


(4.48) 


In  the  dilute  approximation  the  magnetic  chemical  potential 
of  Fe  versus  temperature  is  shown  in  Figure  4.5.  The 
magnetic  chemical  potential  of  Fe  is  not  strongly  dependent 
on  the  concentration  of  solute  X  for  a  typical  impurity 
concentration  ,  such  as  phosphorus  at  0.0001  atomic 
f ract ion , 


Fe 

u 

mag 


GFe  (T)  ,  for  small  XY,  Ta  =  T 
mag '  '  X 


(4.49) 


To  observe  equation  (4.48)  as  a  function  of  magnetic 
parameter,  the  solute  concentration  must  be  high,  as  in 
Figure  4.6.  As  may  be  seen,  the  magnetic  chemical  potential 
of  Fe  increases  with  decreasing  magnetic  parameter . 

Likewise,  the  chemical  potential  of  X  in  alpha  Fe-X  is 


given  by 


•  n  f  n 
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Figure  4.5 

The  magnetic  chemical  potential  of  Fe  in  alpha  Fe- 
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X 

P 


X 

u 

mag 


X 

^  nm 


°g*  +  [i-xy]2iFeX 

nm  u  XJ  nm 


+  RT£nX 


X 

Mmag 


[l-M][G^g(T“)  -  ATXX/S^g] 
+  [l-Xx-MXx+Mx2]ATx[S^g(T“) 


uSFe  j 
mag J 


(4.50) 


(4.51) 

The  magnetic  chemical  potential  of  X  versus  temperature  is 
shown  in  Figure  4.7  for  M  =  1  and  varying  magnetic 
parameter.  Again,  decreasing  magnetic  parameter  increases 
the  potential  curve  and,  for  a  dilute  solution, 


pX  „  =  AT  y [ S  Fe  (T)  -  uSFe  ] 
mag  XL  magv  '  magJ 


(4.52) 


Likewise,  Figure  4.8  shows  the  magnetic  chemical  potential 
of  X  for  varying  solute  concentration. 

The  behavior  of  the  chemical  potential  curves  is  as 
expected.  A  small  addition  of  solute  X  to  alpha  iron  makes 
little  difference  to  the  average  iron  atom  and  therefore  has 
little  effect  on  the  chemical  potential  of  Fe.  However,  each 
solute  X  atom  disturbs  the  alpha  iron  lattice  locally  and 
there  is  a  tendency  for  rejection  of  X  from  the  lattice. 

This  increases  the  chemical  potential  of  X. 

Nishizawa  et  al.  [44],  have  conjectured  that  M  =  0  may 
be  used  for  magnetic  solute  such  as  cobalt  and  M  =  1  may  be 
used  for  non-magnet ic  solute  such  as  phosphorus.  However,  M 
may  be  temperature  and  concentration  dependent.  For  M  =  1, 
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Figure  4.7 

The  effect  of  the  magnetic  parameter  on  the  chemical  potential  of 
(A  =  -2000  K,  B  =  -1000  K,  C  =  0  K.  D  =  1000  K). 
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the  chemical  potentials  are  summarized  below: 


1 . 

Fe 
y  nm 

orFe 

nm 

• — i 

C\J  X 
X 

4- 

Fp  y 

„„  +  RTMl-Xy) 

nm  X 

2. 

Fe 

u 

=  G^6 

(Ta) 

-  ATVXV  uSFe 

mag 

mag 

X  X  mag 

-  Xv[l-XY]ATv[SFe  (Ta)  -  USFe  ] 
XL  XJ  XL  magv  '  magJ 


3. 

X  _  or X  ,  r.  Y  -1 2 T  FeX 
ynm  -  Gnm  +  [1‘XX]  'nm  '  + 

RT£nXv 

X 

4. 

^nm  ’  ^-xx^VSmag<Ta> 

-  uSFe  ] 
mag 

(4.53) 

4.4  Limit  of  Solubility  for  X  in  alpha  Fe-X  Solid  Solution 

A  change  in  the  chemical  potentials  due  to  the  onset  of 
ferromagnetism  must  affect  the  limit  of  solubility  of  a 
solute  in  alpha  iron.  The  general  case  is  considered  first. 


4.4.1  Alpha  Fe-X  Solid  Solutions 

The  limit  of  solubility  for  a  solute  X  in  alpha  Fe-X 
occurs  when  a  second  phase,  FeaX^  (denoted  0),  is  in 
equilibrium  with  the  alpha  Fe-X  solid  solution.  The  a  and  b 
are  defined  as  fractions  such  that  a  +  b  =  1. 

Figure  4.9  shows  a  schematic  of  free  energy  versus 
concentration  for  two  phases  in  equilibrium  at  a  particular 
temperature.  The  dashed  line  represents  a  'tie'  line  that 
determines  the  concentration  of  solute  in  each  phase.  If  X^ 
is  the  concentration  of  X  in  the  alpha  phase,  then 
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Figure  4.9  A  schematic  of  free  energy  for  two  phases  in 
equi librium. 


76 


G“  -  VX  +  [l-xx]uFe  ,4.54) 

and  when-  equations  (4.34)  and  (4.35)  give  the  free  energy  of 
the  solution,  the  '  tie'  line  is  given  by 

r  r  X  Fe Y  F e 

«  “  Lu  -  M  jxv  +  u 

X  (4.55) 

Therefore,  assuming  the  second  phase  is  stoichiometric,  i.e. 

the  atomic  fraction  of  X  in  the  second  phase  is  b, 


G 


0 


/  X  Fe \ .  Fe 
(m  -  li  )b  +  y 


(4.56) 


and  the  condition  for  equilibrium  of  1  mole  of  9  is  given  by 


ap 


Fe 


+ 


(4.57) 


where  G  is  the  molar  free  energy  of  the  second  phase. 

The  following  formalism  was  developed  by  Ko  and 
Nishizawa  [29]  for  alpha  iron.9  For  a  dilute  solid  solution, 
where  X  <<  1 ,  a  dilute  approximation  may  be  applied, 

A 


Fe 


°GFe  + 
nm 


GFe 

mag 


(T) 


X 

u 


°GX  +  IFeX  +  RT£rtXY  +  AT  [SFe  ( T ) - US Fe  ] 
nm  nm  X  X  mag  magJ 


and  equation  (4.57)  may  be  written  as 


(4.58) 


9 Recent  1 y ,  solubility  in  alpha  cobalt  has  also  been 
thermodynamically  analysed  [48]. 
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g9  "  +  «E:9CTI3 

+  b{°GnXm  +  1  nm'  +  RT^*x  +  (T) 


+  ATY[SFe  (T)  -  USFe  ]} 
X  mag  magJ 


(4.59) 


Rearranging  this  equation  gives  the  limit  of  solid 
solubi  1  i-t y ,  ,  as 

A 

XG  =  exp(-[Q  +Q  ] / b  RT ) 

X  Lynm  ymagJ  ' 

Q  =  a°GFe  +  b[°GX  +  IF®X]  -  G9  =  C  +  DT 

ynm  nm  nm  nm 

^mag  =  ^ 1 ' bM]Gm!g < T >  +  bATxtSmag(T)  '  “s«g3  (4.60) 

The  behavior  of  the  magnetic  free  energy  of  solution,  Q 
versus  temperature,  is  shown  in  Figure  4.10  for  M  =  1, 
b  =  0.5  and  varying  magnetic  parameter.  The  magnetic  free 
energy  of  solution  is  approximately  linear  with  temperature 
above  the  Curie  temperature  as  would  be  expected  for  a 
non-magnetic  situation.  However,  this  is  not  true  near  and 
below  the  Curie  temperature;  there  is  a  wide  variation  of 
the  magnetic  free  energy  of  solution  with  magnetic 
parameter.  A  decrease  in  the  magnetic  parameter  increases 
the  magnetic  free  energy  of  solution.  As  the  free  energy  of 
solution  increases,  the  limit  of  solubility  for  the  solute 
will  have  to  decrease. 

Table  4.2  shows  calculated  magnetic  parameters  for  various 
elements  when  added  to  Fe  [44].  Figure  4.11  shows  the 
corresponding  Curie  temperature  for  selected  alpha  Fe-X 
bi nar i es . 


. 
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Table  4.2  Magnetic  parameters  of  solute  in  alpha  iron 
calculated  from  solubility  data  [44] 


Element 

Magnetic  Parameter 

(K) 

Be 

-1800 

Mn 

-1300 

P 

-1200 

Cu 

-1000 

Cr 

-800 

Zn 

-630 

Ni 

-400 

Mo 

-300 

Ti ,Sn,Sb,W 

0 

Co 

+  1020 
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F igune  4.10 

The  effect  of  the  magnetic  parameter  on  the  free  energy  of  so  lution  (A  =  -2000 
B  =  -1000  K,  C  =  0  K,  D  =  1000  K) . 
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The  behavior  of  /n(X^)  versus  inverse  temperature, 
related  to  Figure  4.10,  is  shown  in  Figure  4.12.  As 
expected,  the  solubility  decreases  with  increasing  magnetic 
parameter,  especially  below  the  Curie  temperature  of  the 
binary  alloy.  That  is,  solutes  that  strongly  decrease  the 
Curie  temperature  (large  negative  magnetic  parameter)  also 
strongly  decrease  the  limit  of  solubility.  The  dashed  line 
represents  the  paramagnetic  solubility  line  for  all  magnetic 
parameters.  This  is  a  situation  where  magnetic  free  energy 
is  linear  with  temperature  and  magnetic  entropy  is  constant. 
Therefore,  the  paramagnetic  behavior  is  the  same  as  the 
traditional  ' non-magnet ic'  formulation.  Solubi 1 i ty  data 
above  the  Curie  temperature  approximates  the  paramagnetic 
limit  of  solubility  but  it  must  be  remembered  that  short 
range  magnetic  interactions  still  exist  and  cause  some 
deviation  from  paramagnetic  behavior. 

The  experimental  determination  of  solid  solubility  at 
low  temperatures  (below  800  C)  requires  a  sophisticated 
investigation  technique  with  a  small  spatial  resolution. 
Furthermore,  to  approach  a  state  of  equilibrium,  very  long 
ageing  at  temperature  is  needed.  However  the  advent  of 
quantitative  Electron  Microprobe  Analysis  ( EPMA )  has 
resulted  in  a  substantial  amount  of  experimental  data  over 
the  last  decade.  Some  of  this  data  is  shown  in  Figure  4.13. 
As  is  evident,  the  solubilities  of  many  solutes  in  alpha 
iron  follow  the  form  shown  in  Figure  4.12. 
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F igure  4.12 

The  effect  of  magnetic  parameter  on  predicted  solid  solubility  in  Fe  (A  =  -2000 
B  =  -1000  K,  C  =  paramagnetic,  D  =  0  K,  E  =  1000  K). 
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Solubility  of  Alloying  Elements  (at  %) 
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F igure  4.13 
Examples  of  solid 


solubi 1 i ty 


for 


X 


in  alpha  Fe. 
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4.4.2  The  Limit  of  Solubility  For  P  in  alpha  Fe-P 

The  experimental  results  for  the  solubility  of 
phosphorus  in  alpha  iron  from  the  present  work  are  shown  in 
Table  4.3. 

These  results  were  obtained  from  concentration  profiles  near 
second  phase  precipitates.  The  growth  of  an  iron  phosphide 
precipitate  requires  a  concentration  gradient  from  the 
surface  of  the  particle  into  the  bulk.  If  it  is  assumed  that 
the  region  immediately  near  the  particle  surface  is  in 
quas i -equi 1 i br i urn,  then  the  solute  concentration  of  the  bulk 
in  this  region  is  the  real  limit  of  solubility  for  that 
temperature.  This  is  especially  true  at  lower  temperatures 
as  the  time  to  reach  true  equilibrium  may  be  prohibitively 
long  and  the  concentration  gradient  becomes  very  steep  near 
the  precipitate  surface.  An  example  of  such  a  concentration 
profile  is  shown  in  Figure  4.14. 

The  experimental  limit  of  solubility  for  P  in  alpha  Fe 
(denoted  by  □)  is  shown  in  an  Arrhenius  plot  in  Figure  4.15. 
Also  shown  is  the  calculated  limit  of  solubility  (solid 
line)  and  the  paramagnetic  limit  of  solubility  (dashed 
line).  This  latter  has  been  approximated  from  the  data  of 
Kaneko  et  al.  [23]  who  proposed  the  following  solubility 
equation  (in  at% ) 

Xp  -  exp( 4 . 97  -  4490/T ) 
v  (4.61) 

This  gives  a  non-magnetic  free  energy  of  solution  of 
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Table  4 


3  Experimental  results  of  bulk  phosphorus  content 


Alloy 

phosphorus 

( a  t  % ) 

T800C 

1.97 

T750C 

1.50 

T700C 

0.95 

T650C 

0.72 

T600C 

0.56 

T550C 

0.43 

T500C  . 

no  result 

Master  1 

32.60 

Master  2 

24.62 

ALPHA  IRON-PHOSPHORUS  .  PHOSPHIDE 
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F igure  4. 14 

An  example  concentration  profile  from  T550. 
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Qnm  "  +  6.0T  J/mole 

From  data  obtained  using  EPMA,  Nishizawa  et  al.  [44] 
proposed 


Q  =  9500  +  5.9T  J/mole 
nm 


(4.63) 


and  the  values  b  =  0.25  and  ATp  =  -1200  K.  The  present 
investigation  gave  slightly  lower  concentrations  than  that 
of  Nishizawa  et  al.  This  may  be  due  to  the  fact  that 
Nishizawa  et  al.  used  a  proportional  constant  to  derive 
concentration  from  the  raw  data  of  X-ray  counts.  This  method 
is  not  as  accurate  as  making  iterative  ZAF  corrections  of 
the  raw  data  as  performed  by  the  EDATA2  program. 

The  calculated  solubility  in  Figure  4.15  makes  use  of 
-1200  K  for  the  magnetic  parameter.  The  experimental  data 
fits  well  below  the  Curie  temperature  but  approaches  the 
paramagnetic  limit  too  quickly.  There  are  several  possible 
reasons  for  this  slight  discrepancy.  The  magnetic  parameter 
was  determined  by  a  least  squares  fit  of 


RT£nXp  +  [IzbM]G^g(T)  =  -Q 

5  b  ~ 


-  ATxEsJ^m-^g] 


(4.64) 


over  the  entire  solubility  region.  The  shape  of  the  magnetic 
heat  capacity  curve  (which  determines  the  magnetic  entropy 
and  enthalpy  curves)  is  different  above  and  below  the  Curie 
temperature.  It  may  be  that  the  heat  capacity  curve 
decreases  much  more  rapidly  above  the  Curie  temperature. 
Alternatively  the  magnetic  parameter  may  be  temperature  and 
concentration  dependent  (a  non-linear  variation  of  Curie 
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temperature).  Furthermore,  very  little  is  Known  about  the 
behavior  of  the  solute  intensity  factor,  M. 

In  general,  this  approach  using  empirical  excess 
magnetic  free  energies  seems  justified.  The  thermodynamic 
formalism  fits  the  experimental  data  quite  well.  The 
magnetic  thermodynamic  functions  are  then  valid 
approximations  of  reality  and  it  is  now  possible  to  apply 
this  information  to  other  properties  of  iron-based  alloys. 


4.5  Vacancy  Concentrations  in  alpha  Iron 

Previous  developments  by  the  author  and  others  as 
described  above  are  now  extended  to  the  case  of  vacancies  in 
alpha  iron. 

4.5.1  A  Hypothetical  Iron  Solution 

Consider  an  alpha  iron  crystal  at  absolute  zero.  This 
solid  state  contains  no  vacancies.  Now  assume  there  exists  a 
hypothetical  alpha  iron,  denoted  A,  with  no  vacancies  at  any 
temperature.  Let  the  free  energy  of  this  material  be 

pA  p A  +  qA 

nm  mag  (4.65) 

Assume  that  vacancies  may  be  treated  like  solute  atoms. 
Therefore,  one  has  an  alpha  Fe-v  binary  solid  solution.  By 
analogy  with  an  alpha  Fe-X  regular  solid  solution,  the  free 
energy  of  alpha  Fe-v  may  be  written  as 
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where 


T3  =  T  -  AT  X 
v  v 


(4.66) 


It  is  well  Known  there  are  relatively  few  vacancies  present 
at  lower  temperatures ;  therefore,in  dilute  approximation,  Xv 
is  very  small  and  the  free  energy  of  pure  iron  approximates 
the  free  energy  of  the  hypothetical  iron.  At  a  temperature 
T,  there  is  a  vacancy  content  such  that  the  system  free 
energy  is  a  minimum.  This  occurs  when  the  derivative  of  free 
energy  with  respect  to  the  vacancy  concentration  is  zero.  In 
the  dilute  approximation,  vacancy  concentration  is  small  and 


3G 


Fe 


°GFe  +  RT£rcX 


V 


nm 


where 


°GFe  =  °GA  +  IAV  .  °GV 
nm  nm  nm  *  nm 


(4.67) 


This  gives  the  equilibrium  vacancy  concentration  as 


Q 


X 


mag 


v 


(4.68) 


The  energies  of  vacancy  formation  are  normal ly  measured  at 
high  temperatures  and  therefore  may  be  considered  as  the 
paramagnetic  energies  of  vacancy  formation,  i.e.  when 
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SFe  (T)  =  USFe 
mag  '  mag 

r Fe  / -r \  u  ~  F  6  -r 

mag-  mag  (4.69) 

The  equilibrium  concentration  of  vacancies  in  alpha 
iron  varies  with  temperature  in  a  manner  similar  to  the 
variation  of  the  limit  of  solubility  for  a  solute.  It  is 
expected  that  the  magnetic  parameter  for  vacancies  is  large 
and  negative,  since  the  formation  of  a  vacancy  distorts  the 
iron  lattice  and  there  are  no  electrons  available  for 
coupl i ng . 


4.5.2  Predicted  Vacancy  Concentrations  from  alpha  Iron 
Sel f-Di f fusion  Data 

It  is  well  Known  that  sel f -di f fus ion  in  a  pure  element 
depends  strongly  on  the  vacancy  concentration  present.10  If 
the  vacancy  concentration  is  affected  by  ferromagnet i sm, 
then  sel f-di f fusion  must  be  affected.  The  anomaly  should 
show  up  in  the  sel f-di f fusion  coefficient.  In  fact  the 
se 1 f -d i f f us  ion  coefficient  for  alpha  iron  has  been 
experimentally  determined  by  various  workers  [53-56]. 

Figure  4.16  shows  an  Arrhenius  plot  of  the  self-diffusion 
coefficient  for  alpha  iron  taken  from  Hettich  et  al.  [56]. 

The  similarity  between  this  figure  and  the  solubility 
limits  in  Figure  4.12  is  striking.  Indeed,  experimental 

10The  motion  of  an  atom  is  equivalent  to  the  motion  of  a 
vacancy  in  substitutional  solutions. 


T  emperature, 
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Figure  4.16 

Experimental  se 1 f-di f fusion  data  for  alpha  Fe  [ 53-56 , respect i ve ly] . 
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values  of  the  sel f-di f fusion  coefficient  for  alpha  iron  may 
be  used  to  determine  a  magnetic  parameter  aTv  for  vacancies 
in  alpha  iron. 

Consider  the  paramagnetic  sel f-di f fusion  coefficient 
for  alpha  iron,  D^,  as  measured  by  Hettich  et  al.  [56] 

Dpe  =  Dpexp(-Qp/RT) 

D  =  0.98  cm^/sec 
q  =  236.400  J/mole 

P  '  (4.70) 

Similarly,  a  paramagnetic  diffusion  coefficient  for 
vacancies  in  alpha  iron  has  been  determined  [57] 


D 


v 


0.5  exp ( -88770/ RT ) 


/sec 


(4.71 ) 


The  sel f-di f fusion  coefficient  may  be  given  [57]  in  terms  of 
the  diffusion  coefficient  of  vacancies  and  the  free  energy 
of  formation  of  vacancies,  AGp 


DFe  =  KDvexp ( - AGp/RT ) 


(4.72) 


where  K  =  0.72,  a  correlation  factor.  This  allows  the 
paramagnetic  free  energy  of  vacancy  formation  to  be 
determined  as  147600  -  8.33T  J/mole.  The  non-magnetic 
version  would  be  147600  -  17.33T  J/mole.  Note  that  equation 
(4.71)  also  shows  that  the  free  energy  of  vacancy  migration 
in  alpha  iron  is  approximately  0.4  times  the  free  energy  of 
formation  (i.e.  the  activation  free  energy  for  vacancy 
diffusion  is  approximately  1.4  times  the  free  energy  of 
formation ) . 


. 
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Therefore  a  magnetic  parameter  for  vacancies  was 
determined  by  a  least  squares  fit  similar  to  that  carried 
out  for  equation  (4.64),  and,  Keeping  in  mind  that  this 
value  is  1.4  times  too  large,  the  magnetic  parameter  for 
vacancies  in  alpha  iron  is  found  to  be  approximately 
-2400  K,  a  large  negative  value  as  expected.  Figure  4.17 
shows  the  magnetic  free  energy  of  vacancy  formation.  Vacancy 
formation  becomes  much  more  difficult  in  the  f er romagnet i c 
region.  Figure  4.18  is  the  corresponding  equilibrium  vacancy 
concentration  curve. 

Therefore  the  individual  vacancy  affects  the  magnetic 
properties  of  iron  as  significantly  as  any  solute.  However, 
there  are  generally  few  vacancies  present  so  that  the  total 
effect  i s  sma 1 1 . 


4.6  The  Thermodynamics  of  Equilibrium  Grain  Boundary 
Segregat ion 

It  is  first  necessary  to  understand  the  thermodynamics 
of  interfaces  (in  this  case  grain  boundaries)  and  of 
equilibrium  segregation.  Guttmann's  thermodynamic  treatment 
of  mu  1 1 i -component  segregation  [5,17,58]  is  based  on  the 
formalism  initiated  by  Defay  and  Prigogine  [59].  The  basic 
hypothesis  is  that  an  interface,  4>,  can  be  considered  as  a 
two  dimensional  (2~d)  phase  in  which  all  the  thermodynami c 
functions  can  be  defined  as  in  the  bulk,  a  •  The  fol lowing 
define  the  necessary  variables  and  functions: 
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F igure  4.17 

Predicted  magnetic  free  energy  of  vacancy  formation  in  alpha  Fe. 
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N|  =  number  of  i  atoms  in  <j) 

N?  =  number  of  i  atoms  in  alpha  (bulk) 
a  =  surface  tension  of  <J> 
ft  =  area  of  4> 

=  G  (  P ,  T  ,a ,  N  t,  N?  )  =  Gibbs  free  energy  at  <(> 

Ga  =  G  (P,T,N“) 

=  F  ( V ,  T  ,cr,  N  N?  )  =  Helmoltz  free  energy  at  <|> 
F a  =  F  (V.T.N?) 

The  Helmoltz  and  Gibbs  free  energies  are  related  by: 


G*  =  F*  -  aft 
Ga  =  Fa  +  PVa 


(4.73) 


The  only  difference  between  a  2-d  and  a  3-d  phase  is 
the  surface  tension  of  the  2-d  phase.  The  surface  tension  of 
<f>  is  associated  with  essentially  zero  thickness  and  the 
extensivity  of  4>  is  described  only  by  the  boundary  area,  ft. 
This  is  equivalent  to  ignoring  the  volume  of  <J>  and  thus  the 
PV  terms  in  the  free  energy  of  <)>.  Also,  the  <J>  is  assumed 
incompressible  in  the  thickness  direction,  the  mechanical 
contributions  to  G^  arising  only  from  the  2-d  term  +aft.  The 
free  energy  of  <j)  is  a  function  of  the  intensive  variable 
which  means  that  the  entropy  at  4>  and  the  chemical  potential 
of  i  at  <j)  have  to  be  calculated  at  constant  surface  tension. 

A  simple  form  of  the  boundary  chemical  potentials  may 
be  assumed  by  analogy  with  bulk  functions: 
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ual  =  %al  (P,T)  +  RUna“ 
p41  =  (P,T,o)  +  RT£«a? 


(4.74) 


Also  by  analogy 

F4  =  IN?  S41 

4  1 


(4.75) 


where  £  is  the  Helmholtz  equivalent  of  chemical  potential 
for  i.  The  partial  molar  area  of  i  is  defined  as: 


i 


(9$7/3N^)p  j 

1  r *  1  ,a,r jii 


(4.76) 


and,  since  the  total  boundary  area  obeys  Euler's  theorem  for 
extensive  quantities: 


ft  =  I  N^w-j 
i 

Likewise,  Gibbs  free  energy  may  be  given  by: 


(4.77) 


,4  .  z  N4w4i 


1 

=  -  aft 

which  allows  u  ^  =  S  ^  G w -j  (4.78) 

Consider  a  system  with  grain  boundaries,  4>.  The  equilibrium 
segregation  and  surface  tension  are  given  by: 


X*  =  f(P,T,X®,X®, - ,X®) 

o  =  g(P,T,x“.x“,....,x“)  (4.79) 

The  conditions  for  equilibrium  between  a  two-  and  a  three- 
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dimensional  phase,  for  all  i,  are 


P 


4>i 


(4.80) 


For  an  ideal  solution  the  activities  are  the  concentrations, 
and 


P 


4>i 


or$i 


au.  +  RTlnX?  =  °yal  +  RUnX? 


(4.81 ) 


However,  for  non-ideal  solutions,  partial  excess  free 
energies,  may  be  introduced  so  that: 


for 


V1'  +  ES*  -  eu, 

i  =  1 , 2 , . . . , n- 1 . 


+  RUnX4  =  %ai  +  EG?  +  RUnX“  (4  ^ 

The  solvent  is  given  by  i  =  n,  so  that: 


o  (in  ,  .  dtp  —  o  an  ^  n  t  & 

'  +  Gn  +  RU,lXh  -  u  +  Gn  +  RU,lXn  (4.83) 

Multiplying  equation  (4.83)  by  «./w  and  subtracting  it  from 

equation  (4.82)  to  eliminate  a  gives: 

w./w 


RT£n( [Xa/X4j  1  n  xf/xf)  =  AG1 
n  n  ii 

4[V" 


-  U). 


o  an  ,  -r 

y  +  b 


n 


0). 


n 


EGa] 
n J 


-  [°i;4'1'  -  0wal  +  EG4  -  E§“]  (4.84) 

Where  AGi  is  the  segregation  binding  energy  for  atoms  i  at 


the  boundary. 

Guttmann  and  McLean  [58]  have  shown  that  for  the  case 


of  regular  behavior  with  solute  site  competition11,  the 
segregation  is  given  by: 
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X?  =  X*  exp ( AG1 / RT ) 

rTH  T“ 

1  +  E  [exp ( AGJ /RT)  -  1]X? 

j=l  J  (4.85) 

For  a  non-competitive  case,  each  solute  species  takes  up  a 

different  set  of  sites  in  the  boundary.  If  Y^1  is  the 

fraction  of  sites  filled  by  solute  i,  then,  in  dilute  bulk 

approximation , 

Y?  =  Y“  exp(AG’/RT) 

1  +  Y“  exptAG’/RT)  (4.86) 

Any  or  all  of  the  solute  species  present  may  saturate 

( Y $  =  1)  the  boundary  under  the  appropriate  conditions.  The 
i 

segregation  behavior  of  any  solute  is  independent  of  the 
segregation  of  any  other  solute  except  through  the  free 
energy  of  segregation.  Examples  of  segregation  curves  are 
given  in  Figure  4.19.  These  curves  were  generated  using  the 
computer  program  developed  in  the  present  work  without 
magnetic  effects  (setting  magnetic  terms  to  zero).  It  should 
be  noted  that  in  a  non-competitive  ternary  system,  the 
segregation  curves  of  both  solutes  will  be  similar  to  the 
binary  behavior  described  in  section  4.7.  It  has  been 
suggested  by  Guttmann  [60]  that  similar  atoms  are 
competitive  (e.g.  metalloid  atoms  such  as  P,As)  and  unlike 
atoms  are  non-competitive.  However,  the  evidence  is  not 

1 1 A 1 1  atoms  at  the  boundary  compete  for  one  set  of  boundary 
si tes . 
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conclusive. 


4.7  Grain  Boundary  Segregation  in  Binary  alpha  Fe-X  Solid 
Solutions 

Guttmann's  analysis  is  now  applied  to  binary  alpha  Fe-X 
solid  solutions  using  the  chemical  potentials  obtained  in 
section  4.3  and  found  in  equations  (4.44)  to  (4.48).  For 
convenience  of  presentation  the  following  simplifications 
will  now  be  applied  (excess  energy  terms  are  assumed  to  be 
non-magnet i c ) 


o  4>i  _  o  <{>i  Es<j> 

^  n  m  ^  13  i 

(4.87) 

4.7.1  Case  I:  Grain  Boundary  Free  Energy  Without  Magnetic 
Terms 

As  a  first  approximation,  the  grain  boundaries  will  be 
considered  non-magnet i c .  That  is,  the  atoms  in  the 
boundaries  cannot  take  part  in  any  magnetic  interactions  and 
there  are  no  magnetic  terms  possible  when  describing  the 
boundary  free  energies.  As  the  magnetic  free  energy  terms 
would  be  directly  dependent  on  the  boundary  solute 
concentrations,  this  simplified  case  allows  direct  analytic 
solutions  which,  despite  their  limitations,  do  indicate 


trends . 
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To  obtain  an  equilibrium  segregation  equation,  the 
following  conditions  are  applied  for  some  temperature,  T: 


<J>X  aX 
y  =  y 

d)Fe  aFe 
y  =  y 


(4.88) 

a 


Concentrations  of  solute  X  in  the  bulk  and  boundary  are 
and  Xvx  respectively,  and,  as  with  equations  (4.82)  and 
(4.83) 


y 


<l>x 


°y(J)X  +  Rlln^  -  acoy 
nm  x  X 


°yaX  +  RT£nX“ 
Mnm  X 


y 


aX 

mag 


=  y 


aX 


(4.89) 


u'fFe  =  °y't'Fe  +  RUk(1-X$)  -  auFp 

nm  Xre 

=  °yaFe  +  RUn(l-X“)  +  y°T® 

Mnm  v  X  mag 

aFe 

=  ^  (4.90) 

The  excess  free  energy  terms  are  taken  to  be  non-magnet i c . 
Multiplying  equation  (4.90)  by  and  subtracting  from 

equation  (4.89)  results  in 

RT£n({[l-X“]/[l-X^]}“x/“FeX^/X“)  -  IAGX  =  AGXm  + 

(4.91 ) 


AG 


X 

nm 


I 


i-o  $Fe 

L  r 


nm 


wFe 


o  aFe 
unm 


<$>X 

nm 


aX 

umag 


,  aFe 
^X  pmag 

wFe 


o  aX1 
unm-* 


(4.92) 
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The  superscript  I  in  the  segregation  free  energies  refers  to 
Case  I.  Superscript  II  will  refer  to  Case  II  (section 
4.7.2)  . 

At  this  point,  it  is  assumed  that  the  atomic  structure 
of  the  grain  boundary  is  not  altered  throughout  the 
segregation  process.  In  other  words,  substitution  is  perfect 
with  both  species  occupying  the  same  area  in  <J> f  that  is: 


“Fe  =  “x 


(4.93) 

Equation  (4.92)  then  becomes  (applying  the  approximation  in 
equation  (4.35)): 


!AGX  =  -MGFe  (T)  +  [  1  -  2MX“] ATyIs!" ®  fTa)  -  ] 

mag  magv  1  X  mag  mag  (4.94) 

The  behavior  of  this  magnetic  free  energy  of  segregation  may 
be  evaluated  by  considering  a  dilute  bulk  approximation, 
where  X  ^  <<  1 ,  and 


Fe 


u.Fe 


l4951 

Figure  4.20  shows  the  behavior  of  equation  (4.94)  with 
temperature,  using  the  previously  established  equations  for 
thermodynamic  functions. 

It  is  irmediately  evident  that  the  magnetic  segregation 
free  energy  has  a  large  value  at  low  temperatures.  With 
increasing  temperature,  there  is  first  a  gradual  decline, 
then  an  increasing  rapid  decline  near  the  bulk  Curie 
temperature  culminating  in  a  minimum  just  above  the  Curie 
temperature.  Above  this  point,  the  curve  approaches  a  linear 


ur  Fe 


105 


Q 

G 

G 

G 

(9 

G 

G 

G 

G 

in 

G 

in 

G 

t'* 

in 

<\i 

<\J 

H 

H 

H 

A3W3N3  33&d  N0I±V03a03S  OVW 


* 


111 

q: 

D 

H 

< 

o: 

uj 

a. 

z 

UJ 

H 


C 

o 


03 

E 

•r— 

x 

O 

£_ 

a 

a 

03 


03 

_c 

cl 

03 

c 


a_ 


C3H0W/r> 


Figure  4.20 

Free  energy  of  segregation  for 
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increase  with  temperature.  This  may  be  seen  from  equation 
(4.94)  by  considering  that 

SFe  (T>>TFe)  =  USFe 
mag v  c  '  mag 

GFe(T>>TFe)  =  -TUSF® 
magv  c  mag 

I Ar X  ,r  vTFe^  _  ,MTucFe 

AG  (  T  >  >  T  ]  *  ^nl  S  / «  qc \ 

mag '  c  mag  (4.9b) 

Increasing  the  bulk  solute  concentration  decreases  or 

slightly  flattens  the  curve  below  the  Curie  temperature. 

In  order  to  appreciate  the  effect  of  the  segregation 
free  energy  on  equilibrium  grain  boundary  segregation, 
consider  the  segregation  equation  (4.84).  For  constant 
temperature,  as  the  free  energy  of  segregation  increases, 
the  boundary  solute  content  also  increases.  This  is 
equivalent  to  displacing  the  segregation  curve  to  the  right 
as  shown  in  Figure  4.21.  To  evaluate  his  equilibrium 
segregation  equations,  Guttmann  set  a  constant  standard  free 
energy  of  segregation.  That  is,  there  is  only  one  type  of 
site  for  segregated  solute  atoms  at  the  grain  boundary. 
Likewise,  in  evaluating  the  magnetic  contribution  to 
segregation  in  the  present  work,  the  non-magnet ic  components 
of  segregation  free  energy  are  set  to  a  constant  value. 

These  assumptions,  that  all  partial  molar  boundary 
areas  are  equal  and  that  the  non-magnetic  segregation  free 
energy  is  a  constant,  are  quite  restrictive.  In  reality, 
there  is  a  'spectrum'  of  available  grain  boundary  sites  for 
solute  segregation.  Each  'type'  of  site  will  lower  (or 
raise)  the  energy  of  the  system  a  particular  amount  when  it 
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Figure  4.21 

The  effect  of  the  segregation  free  energy  on  intergranular  solute  concentration 
(A  =  non -magnet ic ,  B  =  magnetic). 
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is  occupied  by  a  particular  species  of  atom.  This  is  Known 
as  the  'spectrum  of  binding  energies'  approach  [61]. 
Unfortunately,  the  spectrum  of  binding  energies  cannot  be 
theoretically  calculated  and  there  is  a  paucity  of  available 
experimental  data  [62].  Furthermore,  competition  for  grain 
boundary  sites  may  occur  if  the  binding  energies  for 
different  solute  species  are  similar.  The  coverage  of  solute 
will  also  vary  from  grain  boundary  to  grain  boundary  as  the 
sites  available  in  the  boundary  will  change  with  the 
structure  of  the  boundary  [63]  which  changes  with  boundary 
mi sor ientat ion ,  etc. 

As  an  example,  consider  the  alpha  Fe-P  system. 

Figure  4.22  shows  a  grain  boundary  concentration  versus 
temperature  plot  for  an  alpha  Fe  -  0.1  at%  P  solid  solution. 
Three  curves  are  presented  for  non-magnet i c ,  paramagnetic, 
and  ferromagnetic  situations  for  the  same  alloy.  It  is  clear 
that  the  effect  of  the  magnetic  contribution  is  to  increase 
segregation  and  its  temperature  dependence  over  a 
commercially  important  range  of  temperatures 
(700  to  1000  K).  The  effect  of  increasing  bulk  solute 
concentration  is  to  raise  the  segregation  curve  and  shift  it 
to  the  right.  If  the  dilute  bulk  approximation  is  not  made, 
so  that  terms  with  powers  of  must  be  included  in  the 
chemical  potentials,  the  effect  is  to  shift  the  segregation 
curve  to  slightly  lower  temperatures . 

The  magnetic  component  of  segregation  free  energy  as  a 
function  of  temperature  for  several  values  of  magnetic 
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Figure  4.22 

Predicted  grain  boundary  solute  concentrations  for  the  alpha  Fe-P  system  according 
to  case  I . 


110 


parameter  is  shown  in  Figure  4.23.  The  corresponding 
boundary  segregation  curves  are  shown  in  Figure  4.24.  for  an 
arbitrary  but  reasonable  value  of  AGpm  =  30  K J /  mole.  It  is 
clear  that  solutes  which  depress  the  Curie  temperature 
(negative  magnetic  parameters)  shift  the  segregation  curve 
to  higher  temperatures .  Likewise,  solutes  which  increase  the 
Curie  temperature  (positive  magnetic  parameters)  shift  the 
segregation  curve  to  lower  temperatures.  The  shift,  as 
indicated,  may  be  hundreds  of  degrees.  The  grain  boundary 
segregation  curves  for  several  bulk  solute  concentrations  is 
shown  in  Figure  4.25  using  the  appropriate  magnetic 
parameter  for  phosphorus.  This  figure  also  displays  the 
effect  of  the  solubility  of  P  in  alpha  Fe  as  shown  by  the 
solid  segment  of  curve  D.  For  the  dashed  lines  the 
solubility  is  ignored.  The  boundary  solute  content  is 
proportional  to  the  bulk  solute  content  and  the  exponential 
of  the  segregation  free  energy  (see  equation  (4.85)).  At 
high  temperatures  the  exponential  term  is  small  which  allows 
solubility  effects  (changing  the  bulk  solute  content  term) 
to  be  apparent.  At  low  temperatures  the  exponential  term 
completely  dominates  the  bulk  term  and  solubility  effects 
are  not  apparent.  Therefore,  changes  in  the  bulk  solute 
concentration  or  the  free  energy  of  segregation  terms  have 
similar  effects  on  the  segregation  curve.  All  these  facts 
are  completely  consistent  with  the  changes  in  the  limit  of 
solid  solubility  as  determined  by  the  magnetic  parameter. 
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Figure  4.23 

The  magnetic  component  of  the  segregation  free  energy  according  to  Case 
(A  =  1000  K,  B  =  0  K,  C  =  -1000  K,  D  =  -2000  K). 
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Figure  4.24 

Predicted  grain  boundary  solute  concentration  according  to  Case  I  (A  =  1000 
B  =  0  K,  C  =  -1000  K,  D  =  -2000  K) . 
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Figure  4.25 

The  effect  of  bulk  solute  content  on  predicted  grain  boundary  solute  concentration 
according  to  Case  1  (A  =  0.00001,  B  =  0.0001,  C  =  0.001,  D  =  0.01). 
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As  is  evident,  Case  I  immediately  points  out  the 
general  trends  of  the  effect  of  ferromagnetism  on 
segregation,  however,  some  inconsistencies  are  present.  The 
aTx.  =  0  K  condition  should  represent  no  magnetic  effect, 
i.e.,  zero  variance  of  the  Curie  temperature,  and  the 
magnetic  contribution  to  the  segregation  free  energy  should 
be  zero.  It  may  be  seen  that  Figure  4.23  shows  a  non-zero 
magnetic  segregation  free  energy  for  the  condition 
aT^  =  0  K.  Also,  above  the  Curie  temperature  the  free  energy 
slowly  increases  for  any  magnetic  parameter,  a  situation 
which  is  clearly  not  reasonable.  Therefore,  Case  II  will  now 
be  considered. 


4.7.2  Case  II:  Segregation  when  Grain  boundaries  may  be 
Ferromagnet ic 

The  assumption  that  the  only  difference  between  a  2-d 
and  a  3-d  phase  is  the  2-d  surface  tension,  a  .  realistically 
demands  that  the  grain  boundary  of  a  ferromagnetic  3~d  phase 
may  itself  be  fer romagnet i c .  Thus,  the  chemical  potentials 
of  the  grain  boundary  species  would  have  magnetic  components 
(which  they  do  not  have  in  equations  (4.87)  and  (4.88)). 
Hence,  one  may  write 


vi 


i  (p  i 


=  u 


nm 


*i 

mag 


o  <J>i 
u 

nm 


+  RUnX] 


4>  . 


00). 


mag 


(4.97) 


where  is  a  magnetic  free  energy  component  of  the 

boundary.  It  is  assumed  that  the  magnetic  chemical  potential 


'  *  '  ■’  ^ 
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in  <J>  may  be  described  as  in  the  bulk,  i.e. 


4>Fe 

u 

mag 


Gmag<T*>  '  "Vx^mag 
-  X*n-MX*]ATx[S^g(T*) 


uSFe  ] 
mag 


(4.98) 


umag  ’  H -M X J g ( T* )  -  MxX*USF*g] 

+  [l-X^-HX*+Mxf]ATx[S^g(T^)  -  US^g]  (4  9g) 

where  T*  =  T  -  aTxxf.  As  in  Case  I,  the  Guttmann  analysis  is 
applied  to  obtain  an  equilibrium  segregation  equation,  i.e. 

for  some  temperature,  T, 

4>X  =  o  *X  +  RUnX*  .  +  *X 

H  Mnm  X  X  mag 

0  X  ,  DT/?  y  0^  i  OtX 
=  u  +  RTcnXw  +  y 
Mnm  X  mag 

aX 

=  P 

.  op,Fe  +  RUn(i.xJ)  -  „Mfe  ♦  y«£| 

-  Xm6  +  RT£«(1-X“)  +  y“^ 

=  yctFe  (4.100) 

After  treatment  as  in  Case  I  and  again  assuming  ^pe=  wx ,  a 
segregation  equation  results: 
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RTMx*[i-x*]/x“[i-xJ])  =  hagx 


=  AGX  +  IJAGX 
nm  mag 

HAGX  =  M[GFe  ( )  -  G Fe  ( Ta )  -  [X^-XY]ATYUSFe  ] 
mag  mag  mag  XXX  mag 

-  [1  -MxJ]AT¥[SFe  (T*)  -  uSFeJ 
X  X  mag  mag 

+  [l-MX“]ATv[SFe  (T“)  -  USFe  ] 

X  X L  mag  mag J  (4.101) 

No  technique  was  found  to  evaluate  this  equation  directly 
but  by  considering  equation  (4.85),  then 


xj  =  x£  exp ( AGX/RT ) 

1  +  x“[exp(AGX/RT)  -  1] 


(4.102) 


is  solved  for  roots  of  x|.  This  may  be  done,  for  example,  by 

the  graphical  bisection  method  or  an  iterative  substitution 

of  X^,  and/or  T.  The  computer  programs  and  techniques 

developed  for  these  numerical  solutions  can  be  found  in 

Appendix  II.  To  better  understand  the  expression  for  the 

magnetic  free  energy  of  segregation  consider  the  two 

bounding  ranges,  high  temperatures  (T  >>  T )  and  low 

temperatures  (T  <<  T  )  .  At  higher  temperatures,  above  the 

Curie  temperature,  the  boundary  and  bulk  solute 

concentrations  are  similar  (X^  ~  X^)  and  since  the  magnetic 

segregation  free  energy  is  the  difference  between  magnetic 

components  of  the  grain  boundary  and  the  bulk  segregation 

free  energies,  (T  >  T^)  will  tend  to  zero  (the 

a  mag  c 

magnetic  effects  are  effectively  zero  above  the  Curie 
temperature ) .  As  the  grain  boundary  becomes  saturated  with 
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solute  (X^  approaches  1)  at  lower  temperatures,  the  shift  in 
Ta  (i.e.  T  -  aT^X^  )  reac^es  a  maximum.  For  typical  values  of 
AT  =  -1000  K,  this  implies 

GFe  ( )  <  GFe  ( Ta ) 
mag v  1  magv 

SFe  ( Ta 1  <  SFe  ( T^)  -  uSFe 
mag  ^  1  magv  '  mag  (4.103) 

Also,  as  X  approaches  saturation,  the  boundary  solute 
concentration  is  much  greater  than  in  the  bulk  (X$  >>  X^, 
thus  the  terms  containing  X*J  and  S  (Ta)  may  be  neglected. 
The  overall  effect  is  that  at  low  temperatures 

“  MtGmag<T*>  '  0T<X,] 

-  ATXUsmag[1+M]  (4.104) 

The  first  term  is  negative  and  proportional  to  US^  > 
however  the  second  term  has  a  much  greater  positive  value. 
This  leads  to  an  increase  in  segregation  free  energy  and 
hence  in  equilibrium  segregation,  the  reverse  being  true  for 
positive  magnetic  parameters. 

Figure  4.26  shows  the  magnetic  contribution  to  the 
segregation  free  energy  as  a  function  of  temperature  for 
various  values  of  magnetic  parameter.  The  non-magnetic 
component  of  segregation  free  energy  was  chosen  as 
40  KJ/mole.  This  value  is  higher  than  that  used  in  Case  I 
because  of  the  constant  positive  (background)  magnetic 
segregation  free  energy  of  approximately  10  KJ/mole  found  in 
Case  I  (see  Figure  4.23).  Again,  a  large  negative  magnetic 
parameter  leads  to  a  larger  segregation  free  energy.  For  the 
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Figure  4.26 

The  magnetic  component  of  segregation  free  energy  according  to  Case  II  (A 
B  =  0  K,  C  =  -1000  K,  D  =  -2000  K) . 


aTx  =  0  K  condition,  the  magnetic  free  energy  is  indeed  zero 
for  all  temperatures  and,  above  the  bulk  Curie  temperature, 
all  magnetic  free  energies  tend  to  zero.  Case  II,  then,  is 
more  realistic  than  Case  I.  The  corresponding  amounts  of 
grain  boundary  segregation  for  Figure  4.26  are  shown  in 
Figure  4.27.  It  is  clear  that  decreasing  values  of  magnetic 
parameter  cause  the  segregation  curves  to  shift  to  higher 
temperatures  and  they  can  also  develop  a  sharp  step  below 
the  bulk  Curie  temperature.  Similar  steps  appear  in  the 
segregation  curves  of  ternary  alloys  [58]  where  increases  in 
segregation  free  energy  are  caused  by  solute-solute 
interaction.  For  magnetic  parameters  large  and  negative,  the 
magnetic  free  energy  of  segregation  curves  may  be  divided 
into  3  regions: 

1.  at  temperatures  above  the  Curie  temperature,  the  free 
energy  is  very  small; 

2.  at  some  temperature  usually  below  the  Curie  temperature 
there  is  a  rapid  increase  in  the  free  energy;  and 

3.  at  lower  temperatures  the  free  energy  increases  with 
decreasing  temperature  in  a  roughly  linear  manner. 

The  effect  of  bulk  solute  concentration  on  the  magnetic 
segregation  free  energy  is  shown  in  Figure  4.28,  and  the 
corresponding  grain  boundary  concentration  curves  in 
Figure  4.29.  Interestingly,  the  shape  of  the  magnetic 
segregation  free  energy  curves  is  reminiscent  of  a 
magnetization  curve,  having  large  values  at  low  temperatures 
and  dropping  rapidly  to  very  small  values  near  the 
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Figure  4.28 

The  effect  of  bulk  solute  content  on  the  magnetic  segregation  free  energy  according 
to  Case  II  (A  =  0.00001,  B  =  0.0001,  C  =  0.001,  D  =  0.01). 
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Figure  4.29 

The  effect  of  bulk  solute  content  on  predicted  grain  boundary  solute  concentrations 
according  to  Case  II  (A  =  0.00001,  B  =  0.0001,  C  =  0.001,  D  =  0.01). 
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order -di sorder  transition  temperature  (Curie  temperature). 
The  effect  of  decreasing  bulk  solute  is  to  'chop  off'  a 
right  hand  portion  of  the  magnetic  segregation  free  energy 
curve  by  introducing  steps  at  decreasing  temperatures,  thus, 
shifting  the  region  of  rapid  segregation  increase  to  lower 
temperatures.  Again,  consider  the  segregation  equation 
(4.85).  It  is  obvious  that  the  changes  in  the  bulk  solute 
concentration  and  the  exponential  term  should  have  similar 
effects . 

The  effect  of  non-magnetic  segregation  free  energy  on 
the  magnetic  segregation  free  energy  is  shown  in 
Figure  4.30.  and  the  corresponding  segregation  curves  in 
Figure  4.31.  These  figures  are  very  similar  to  Figures  4.28 
and  4.29  for  magnetic  segregation  free  energy.  • 

The  present  analysis  predicts  that  small  amounts  of 
solute  with  large  negative  magnetic  parameters  have  a  large 
effect  on  the  equilibrium  segregation  behavior.  Thus  grain 
boundary  segregation  behavior  in  iron-based  alloys  is 
different  than  segregation  in  non-magnetic  alloys.  In 
general,  as  many  important  impurities  in  iron-based  alloys 
are  non-magnetic  with  large  negative  magnetic  parameters 
(see  Table  4.2),  the  temperature  range  over  which  the  grain 
boundaries  are  saturated  extends  to  higher  temperatures  than 
would  be  expected  from  purely  non-magnetic  considerations. 

At  these  higher  temperatures,  impurity  diffusion  is  easier 
and  equilibrium  segregation  is  achieved  more  rapidly.  As 
intergranular  embrittlement  of  many  materials  has  been 
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Figure  4.30 

The  effect  of  the  non-magnetic  segregation  free  energy  on  the  magnetic  segregation 
free  energy  according  to  Case  II  (A  =  20  KJ/mole,  B  =  30  KJ/mole,  C  =  40  Kd/molef 
D  =  50  KJ/mole) . 
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Figure  4.31 

The  effect  of  the  non-magnetic  segregation  free  energy  on  the  grain  boundary  solute 
concentration  according  to  Case  II  (A  =  20  Kd/mole,  B  =  30  Kd/mole,  C  =  40  Kd/mole, 
D  =  50  Kd/mole ) . 


126 


linked  to  grain  boundary  impurity  segregation,  the  present 
analysis  indicates  magnetic  materials  are  likely  to  be  more 
susceptible  at  higher  temperatures  to  embrittlement  than 
equivalent  non-magnetic  materials. 

Approximations  of  Case  I  and  Case  II  for  the  alpha  Fe-P 
system  are  shown  in  Figure  4.32.  This  figure  emphasizes  the 
major  differences  between  Case  I  and  Case  II,  that  is,  Case 
II  may  have  a  step  in  the  curve  and  the  Case  I  curve  is 
displaced  to  higher  temperatures.  The  displacement  is  due  to 
the  fact  that,  for  Case  I,  the  magnetic  segregation  free 
energy  has  an  unrealistic  positive  component  present  as  seen 
in  Figure  4.23.  The  step  in  the  Case  II  curve  is  believed 
due  to  a  magnetic  transition  which  will  now  be  discussed. 


4.8  The  Magnetic  State  of  the  Boundary  Phase 

An  interesting  outcome  of  the  above  analysis  arises 
from  the  solute  dependence  of  the  Curie  temperature. 
Boundaries  rich  in  solute  (with  a  large  negative  magnetic 
parameter)  can  achieve  transition  to  the  paramagnetic  state 
as  the  Curie  temperature  of  the  boundary  falls  below  the 
system  temperature.  Thus  situations  can  arise  such  that  as 
the  temperature  falls,  the  boundary  phase  transforms  from 
paramagnetic  to  ferromagnetic  near  the  bulk  Curie 
temperature  then  at  a  lower  temperature  segregation  causes 
the  boundary  Curie  temperature  to  fall  low  enough  that  the 
boundary  phase  again  becomes  paramagnetic. 
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and  Case  II  segregation  curves  for  the  alpha  Fe-P  system. 
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The  grain  boundary  Curie  temperature  is 


Tc  =  Tce  +  iT“X 


<p 

X"X 


(4.105) 


Fe 

where  T  is  the  Curie  temperature  of  iron  (1043  K).  Thus  at 
a  temperature  T,  the  grain  boundary  will  be  ferromagnetic  if 

1  [T  -  T^e],  ATX  negative 

atZ 


xj  >  [T  -  T^], 

IT 


AT ^  positive 


x  (4.106) 

At  any  temperature  a  maximum  solute  content  Xj*  (or  critical 

solute  content)  exists  beyond  which  the  grain  boundary  phase 
« 

will  be  paramagnetic.  A  line  representing  X^*  versus 
temperature  is  shown  for  the  Fe-P  system  in  Figure  4.32. 
Shown  in  Figure  4.29  is  a  more  general  case  at  several  bulk 
solute  contents.  It  can  be  seen  that  in  this  case,  for 
Xx  =  0.001,  the  boundary  is  always  paramagnetic  whereas  for 
Xx  =  0.0001  the  double  transformat  ion  occurs.  It  is  noted 
that  a  step  in  the  segregation  curve  exists  only  when  the 
Xf  line  crosses  the  segregation  curve,  that  is,  when  a 
second  magnetic  transition  takes  place.  However,  a  second 
magnetic  transition  may  occur  without  a  step  occurring  in 
the  segregation  curve  for  small  negative  magnetic  parameter. 

A* 

The  segregation  curves  seem  to  indicate  that  the  X^  line 
crosses  at  or  near  an  inflection  point  in  the  segregation 
curve.  Increasing  the  bulk  content  has  the  effect  of 
decreasing  the  temperature  range  for  which  the  grain 
boundary  may  be  ferromagnetic. 


* 
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The  information  in  figure  4.28  may  be  displayed  in  a 
different  manner  which  will  be  more  useful  in  subsequent 
sections.  Figure  4.33  shows  equilibrium  grain  boundary  Curie 
temperatures  versus  the  system  temperature.  The  conditions 
for  a  ferromagnetic  boundary  would  be 

T*  >  T 

C  (4.107) 

The  dashed  line  separates  the  ferromagnet ic  and  paramagnetic 
regions . 


4.9  The  Thermodynamics  of  Dilute  alpha  Fe-M-X  Solid 
Solut ions 

The  theories  and  formalisms  previously  applied  to  alpha 
Fe-X  binary  solutions  are  now  used  to  describe  alpha  Fe-M-X 
ternary  solutions.  Although  Nishizawa  et  al.  have 
empirically  calculated  the  free  energy  and  solid  solubility 
specifically  for  Fe-Co-X,  the  author  now  presents  a  general 
formalism  for  alpha  Fe-M-X  systems. 

The  total  Gibbs  free  energy  function,  Ga,  for  an  alpha 
Fe-M-X  solid  solution  will  be,  by  definition, 


G“(T)  =  Ha(T)  -  TSa(T) 


(4.108) 


As  with  binary  solutions,  the  total  free  energy  functions  of 
ternary  solutions  may  also  be  analysed  by  their  component 
parts.  That  is,  it  is  assumed,  by  analogy,  that 
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Figure  4.33 

The  effect  of  bulk  solute  content  on  the  boundary  Curie  temperature. 
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Ga(T )  =  Ga  (T)  +  Ga  (T) 
v  '  nmv  '  magv  ' 


(4.109) 


The  non-magnet ic  component,  GJjm,  may  be  theoretically 


described  by  way  of  the  ordinary  regular  solution 
approximation  (the  quasi -binary  approximation).  The  regular 
solution  approximation  gives  the  following: 


.ex  _  x  OqFg  x  OpM  ^  o  q  X 
'nm  Fe  nm  M  nm  x  nm 

+  X..[1-Xm]IF®M  +  Xy[1-Xx]IF®X  +  X,,X  If™ 
M  M  n  m  X  X  nm  M  X  n m 

+  RT[Xpe£nXFe  +  XM£nXM  +  Xx£nXx] 


(4.110) 


where 

0G^m  is  the  non-magnetic  standard  free  energy  for  i, 

I1J  is  the  non-magnetic  interaction  coefficient  between 
nm 

i  and  j, 

X.|  denotes  the  concentration  of  impurity  i,  and 

^denotes  the  concentration  of  a  transition  metal  M. 

A  descriptive  modification  of  the  magnetic  free  energy 
for  the  ternary  solution  is  based  upon  the  modification  used 
for  the  binary  solution.  As  in  the  binary  analysis,  the 
Curie  temperature  of  a  solid  solution  is  assumed  to  vary 
linearly  with  solute  concentration.  Therefore,  by  dissolving 
atomic  fractions  and  xj?j  of  solute, 


TFe  +  ATyXv 
c  xx 


+  atmxm 


where  aTu  is  the  magnetic  parameter  for  M. 


(4.111) 
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The  magnetic  free  energy  of  alpha  Fe-M-X  may  be 
obtained  by  modifying  the  free  energy  of  alpha  Fe  (or  alpha 
Fe-X)  in  such  a  manner  that  the  alpha  Fe-M-X  magnetic  free 
energy  curve  maintains  the  same  position  with  respect  to  a 
new  Curie  temperature  as  the  alpha  Fe  (alpha  Fe-X)  free 
energy  curve  maintains  to  the  old  Curie  temperature.  The 
displacements  may  be  obtained  in  the  manner  described  by 
equations  4.22  to  4.25,  such  that 

T“  =  T  -  [*Txxx  +  fiTriXM]  (4.112 


Ga 

mag 


(T) 


GFe 

mag 


(T“)  -  [ATxXx+ATMXM]US^g 


(4.113) 


The  influence  of  the  solutes  upon  the  size  of  the  magnetic 
component  [41]  introduces  two  constant  solute  intensity 
factors  Mx,  M^,  so  that 


Wg  *  n-MMXM][l-MxXx][G^g(T“)-[ATxXx+ATMXH]uS-a] 


usFe 

X  "  X  ' 1-1 '  M  "  M J  bmagJ  1  14) 


If  aTq  =  ( aTx  Xx  +  aTmXm),  equation  4.114  is  the  form  for  the 
magnetic  free  energy  of  dilute  ternary  solid  solutions 
implicit  in  the  work  of  Mishizawa  et  a  1 . [ 44 ] .  Species  X  and 
M  are  considered  non-magnet ic  and  are  assumed  not  to 
interact  magnetically.  The  magnetic  interaction  terms  may  be 
defined  as  in  equations  (4.28)  to  (4.33)  to  give  the 
interaction  per  unit  mole  of  solute. 
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[ATV  +  AT..]  [S Fe  (Ta)  -  uSFe  ] 
X  MJLmagv  '  magJ 


(4.115) 


This  gives  the  alpha  Fe-M-X  total  free  energy  as 


G“(T)  =  X 


or  Fe 
Fe  bnm 


oJ\ 


0  J 


+  X„  G  +  Xu  G 


M  nm 


X  nm 


+  +  +  wZ 

+  RT[XFeT«XFe  +  XM£nXM  +  XYtnXy] 


M  M 


-  [l-'-'MXH^1-MXXx][GIag(T“)-[ATXXX+ATMXr^Us;eag] 


u,.Fe 


(4.116) 


4.10  Chemical  Potentials  For  Alpha  Fe-M-X 

The  chemical  potential,  jj  1 ,  of  species  i,  in  alpha 
Fe-M-X  is  described  in  the  normal  manner,  with  magnetic  and 
non-magnet ic  components,  so  that, 


Fe 

y 


x 

y 


M 

y 


Ga 

-  X  y  3  G 

d  X 

Fe 

+  3Ga 

3XX 

Fe 

+  3  Ga 

3Xm 

XjjSG01 


3X 


M 


(4.117) 


As  for  the  binary  case  in  section  4.3  various  terms 
must  be  evaluated  to  obtain  a  useful  chemical  potential: 
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i .  a°G1  =  o 

nm 

ax, 

J 


2. 


3l1J  =  0 
nm 

ax, 


3.  axFe  -  -i 


4. 


5. 


6. 


ax 


i^Fe 


9  X  i  ^  F  e  =  0 

ax . 

J 

3GFe  (Ta)  = 
mag  _ 

axi 


AT  .  SFe  ( T° ) 
i  mag 


M.  (X  .  ,T)  =  constant 


(4.118) 


The  partial  differential  of  the  total  free  energy  of 
alpha  Fe-M-X  with  respect  to  the  bulk  concentration  of  X  is 


rv  n  CL 

a  o 


FeX 


=  -°GFe  +  °GX  +  [1 -2XV]I 
nm  nm  nm  X  nm 


ax 


+  X„IMX  +  RT[£nXv  -  £nXcJ 


M  nm 


X 


Fe 


(4.119) 


3G 


a 


u.Fe 


mag 


ax 


-  Mx[l-MMXH][GF®g(T“)-[ATxXx+ATMXM]uS^g] 
+  tl-HMX„][l-HxXx][Smag(T  )-  smag!laTx 


(4.120) 


The  partial  differential  of  the  total  free  energy  of  alpha 
with  respect  to  bulk  concentration  of  M  is 


135 


9Ga  =  -°GFe  +  °GM  +  [l-2XM]IFeM 
nm  nm  nm  L  MJ  nm 


ax 


+  XvIfJX  +  RT[£nXu  -  lnXrJ 
X  nm  M  Fe 


(4.121) 


3Gmag  ■  -  HMf1-MXXX^Gmag(Ta)-^TxXX+ATMXH]Usm®g] 
M 


ax 


mag 

+  [i-hmxh][i-mxxx][s^„(t01)-us!!Jat 


mag 


mag 


mag 


M 


The  chemical  potential  of  Fe  is 


(4. 122) 


Fg  _  or  Fe  ,  y^^Fg  .  v  2 .  F  6  y  v  t  FI  X  ,  n  ~r  d  v 
y  -  G  +XVI  +  X,.  I  -  XMXVI  +RTcnXr 

nm  nm  X  nm  M  nm  MX  nm  Fe 


(4.123) 


vmag  =  ^ -MMMXXMXX^S^gTfa)-[ATxXx+ATMXH]uS^g] 

-  n-MMXH][l-MxXx][&TxXx+ATMXM][Smag(Ta)-  Smag] 


The  chemical  potential  of  X  is 


(4.124) 


nm 


orX 

nm 


+  [l-x  ]2I^®X  -  X..[l-Xv]IFeM 
X  nm  M  XJ  nm 

MX 


+  XM^  "XX  ^  nm  +  RUnhe 


(4.125) 


[l-HX+MHHXXM-VxXHXX^Gmag(Ta)-tATXXX+ATMXM^Smag] 
+  n-MMXH][l-MxXx][ATx[l-Xx]-ATMxM][SF®a(Ta)-uSF!„] 


ma  g 


mag 


(4.126) 


The  chemical  potential  of  M  is 


136 


M 

^  nm 


0 G ^  -  x  r 1  - X  lT^e^  +  n  y  i2TFeM 

h  “  XyU  XMJInm  +  *-l-XM3  Inm 


nm 


+  XX  ^  ”XM^  1  nm  +  RT^X 


M 


(4. 127) 


M 

‘mag 


Fe 

maq 


[|-Mii*B«"x,x-"»M<xn<x1[G«j(T”)-IaTx,i*aT»<Hl"s 
*  [,-"«,HH'-»IX<][»T,[l.X»].aTIX,][S^  (T«)-“s::„] 


mag 
(4.128) 


For  a  dilute  solution,  the  magnetic  chemical  potential  of  Fe 
versus  temperature  is  essentially  the  same  as  shown  in 
Figure  4.5.  The  magnetic  chemical  potential  of  X  in  a 
ternary  has  essentially  the  same  behavior  with  respect  to 
its  magnetic  parameter  as  in  a  binary  solution  (see  Figure 
4.7).  The  effect  of  the  magnetic  parameter  of  M  depends 
largely  on  the  concentration  of  M  but  as  seen  in 
Figure  4.34,  decreasing  magnetic  parameter  increases  the 
chemical  potential.  The  effect  of  the  concentration  of 
solute  M  on  the  chemical  potential  of  X  is  shown  in 
Figure  4.35.  Increasing  solute  content  decreases  this 
potent i a  1 . 


4.11  The  Solid  solubility  of  X  in  Alpha  Fe-M-X  Solid 
Solut ions 

The  limit  of  solid  solubility  of  solute  X  in  alpha 
Fe-M-X  occurs  when  FeaXb  or  MaXb  phases  are  in  equilibrium 
with  the  alpha  Fe-M-X  solid  solution.  Solute  M  wi 1 1  be 
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considered  a  transition  metal  element  such  as  Mn,  Cr,  Mo, 
and  V.  Solute  X  wi 1 1  be  considered  an  impurity  element  such 
as  P,Sb,  and  As.  When  dealing  with  systems  exhibiting  temper 
embrittlement  it  has  been  noticed  [64]  that  the 
intermetal  1 ic  compounds  FeaXb  and  MaXb  have  the  same  crystal 
structure  and  that  complete  miscibility  exists  between  them. 
Therefore  the  compound  beta  ( Fe^M-] _y ) aXb  (where  a  +  b  =  1) 
is  in  equilibrium  with  the  alpha  Fe-M-X  solid  solution.  As 
shown  by  Hillert  [65]  and  by  Guttmann  [64],  the  conditions 
of  equilibrium  between  the  solid  solution  alpha  and  the 
compound  beta  are  obtained  by  treating  beta  as  a  regular 
solution  of  the  two  binary  compounds  FeaXb  (denoted  01)  and 
M  X  (denoted  02)  with  concentrations  y  and  1-y,  that  is 

3  D 

G91  =  °g91  +  a  RT£ny 
G02  =  °G02  +  aRUn(l-y) 

(4.129) 

At  equilibrium,  the  condition  set  out  in  equation  (4.57) 
must  be  applied  such  that 


where 


Fe 

ay 

M 

ay 

i 

y  = 


+  by 

+  bpx 

o  i 

y 


x 


+ 


-  G61 
G02 

RTinX. 


(4.130) 


Equations  4.129  and  4.130  are  combined  to  give  the  system  of 
equations 
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a°pFe  +  aRUrtXFe  +  b°pX  +  bRT£nXx  =  °G01  +  aRT-Gny 

a°uM  +  aRT£nXM  +  b°pX  +  bRT-GnX..  =  °G02  +  aRTGn(l-y) 

(4.131) 


Ar01  or91  o  Fe  ,o  X 

A  G  =  G  -dp  -bp 

.r92_or02  o  M  ,o  X 

AG  -  G  -dp  -  bp 


(4.132) 


To  solve  these  equations  for  the  solubility,  Xs,  it  is 

X 

necessary  to  eliminate  y  and  assume  a  dilute  solution.  The 
exponentials  of  equations  (4.131)  are  taken  and  y  separated 
out  to  give 


y  =  XFeXx/aexp(-AG01/aRT) 
1-y  =  XMX^/aexp(-AG02/aRT) 

Eliminating  y  gives 


(4.133) 


1  =  Xv/a[XF  exp(-AG91/aRT)+X.,exp(-AG02/aRT)]  , 

x  Fe  M  (4.134) 

Applying  the  dilute  bulk  approximation,  where  the 
concentration  of  X  is  small,  gives 


Xx  =  exp(AG01 /bRT) _ 

[l-XM+XMexp(-AG'/aRT)]a/b 

AG'  =  AG02  -  AG01 


(4. 135) 


(4.136) 


* 
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The  numerator  of  equation  (4.135)  denotes  the 
solubility  limit  of  X  in  alpha  Fe  in  the  absence  of  M  and  is 
exactly  equivalent  to  equation  (4.60).  It  is  immediately 
obvious  that  both  aG9^  and  aG'  have  terms  of  excess  free 
energy  that  are  magnetic  in  nature  such  that: 


AG91  =  °Gei  .  [a/e  +  bpX  ]  -  [aMFe  +  buX  ] 

L  nm  pnmJ  L  Mmag  mag 


(4.137) 


AG*  =  °G92 


or91  r  M  Fe  ,  M  Fe 

G  -  alp  -  y  +  y  -  -  y 

L  nm  nm  mag  mag 


(4.138) 


Applying  the  dilute  bulk  regular  solution  approximation,  the 
magnetic  component  of  equation  (4.137)  becomes 


-AG®L  =  [  1 -bMX  +  b  VxXM^Gma g ( T“ 5  ‘ ] 


ur  Fe 


mag 


-  [l-MMxM][bATx-ATM][smag(T  )~  smagJ  (4.139) 
Likewise,  the  magnetic  component  of  equation  (4.138)  is 


ur  Fe 


AG 


mag 


=  -  a  3G 


a 


3  X . 


*  M  (4.140) 

Which  is  equivalent  to  a  constant  times  equation  (4.122)  and 


AGmag  =  aHM[Gmag(T“)  ‘  ATHXMUsmag] 

-  a[1'MMXH^ATM^Smag^T  ^  *  Smag^  (4.141) 
where  Ta  =  T-aT^X^.  Under  these  circumstances  the  limit  of 
solid  solubility  for  X  may  be  obtained  from  equation  (4.135) 
analytical ly . 

The  effect  of  the  magnetic  parameter  of  the  transition 
metal  alloying  element,  aTm,  on  phosphorus  solid  solubility 
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in  alpha  Fe-O.OIM-P  is  shown  in  Figure  4.36.  As  can  be 
noticed,  the  general  shape  of  the  binary  solubility  curve 
(see  Figure  4.12)  is  retained.  However,  a  large  negative 
magnetic  parameter  of  M  further  decreases  the  solubility 
below  the  Curie  temperature  significantly.  An  increase  in 
the  amount  of  transition  metal,  M,  present  also 
significantly  decreases  solubility  of  P  in  alpha  Fe-M-P  as 
is  shown  in  Figure  4.37. 

The  two  previous  figures  are  concerned  with  magnetic 
interactions.  It  is  of  interest  to  also  consider  chemical 
interaction  among  the  elements  as  in  the  work  of  Guttmann 
[58].  He  defined  a  preferential  interaction  coefficient 
between  M  and  X,  r  ,  such  that,  in  the  paramagnetic  limit, 


-AG'  =  abr 


(4.142) 


A  large  positive  I'  denotes  a  strong  attraction  between  M 
and  X  atoms.12  The  effect  of  I'  on  P  solubility  in  alpha 
Fe-O.OIM-P  is  shown  in  Figure  4.38.  As  expected,  large 
positive  I'  drastically  reduces  the  maximum  amount  of  P  in 
solution.  In  other  words,  M  atoms  'scavenge'  P  atoms  and 
form  some  type  of  complex.  The  definition  of  aG'  allows  for 
its  calculation  from  experimental  data.  Guttmann  determined 
I'  for  various  transition  metal-  phosphorus  combinations  as 
shown  in  Table  4.4 

by  making  use  of  the  experimental  solubility  data  of  Kaneko 
[23],  Figure  4.39.  Table  4.4  also  contains  magnetic 
1 21'  wi i i  be  more  fully  explained  in  section  4.13. 
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Table  4.4  Interaction  parameters  for  M-P  calculated  from 
solubi 1 i ty  data  [ 64 ] 


ement 

Interaction 

parameter 

Magnetic  Parameter 

I'  (KJ/mole) 

(K) 

Ni 

30.1 

-400 

Mn 

66 . 9 

-1300 

Cr 

100.3 

-800 

V 

143.8 

1000 

W 

170.7 

- 

Mo 

184.0 

-300 

Ti 

221.0 

- 

Zr 

256.4 

- 
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Figure  4.36 

The  effect  of  the  magnetic  parameter  of  M  on  phosphorus  solid  solubility  in  alpha 
Fe-O.OIM-P  (A  =  -2000  K,  B  =  -1000  K,  C  =  0  K,  D  =  1000  K). 
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Figure  4.38 

The  effect  of  I'  on  the  solubility  of  phosphorus  in  alpha  Fe-O.OIM- 
(A  =  100  Kd/mole,  B  =  75  Kd/mole,  C  =  50  Kd/mole,  D  =  25  Kd/mole). 
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T(°C) 


104/T(°K) 


Figure  4.39  The  experimental  solid  solubility  of  phosphorus 
for  selected  ternaries  [23]. 
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parameters  (where  available)  for  the  metals.  There  is  no 
direct  relationship  between  I'  and  the  magnetic  parameter. 

When  considering  Figure  4.39,  it  should  be  noted  that 
all  solubility  determinations  lie  above  the  Curie 
temperature.  Therefore,  a  calculated  aG'  represents  a  value 
close  to  its  paramagnetic  limit.  However,  there  are  still 
some  short  range  magnetic  interactions  taking  place.  In 
general,  the  magnetic  component  of  aG'  raises  the  total 
value  of  aG' .  The  calculated  values  of  V  should  be 
therefore  considered  high.  A  paramagnetic  I#  being  defined 
constant  is,  in  fact,  a  type  of  zeroth  approximation  with 
respect  to  temperature.  A  future  consideration  may  be  the 
determination  of  an  I'  that  is  linear  with  system 
temperature  (i.e.  a  first  approximation).  In  the  present 
work,  Guttmann's  values  of  I'  will  be  used  when  combining 
magnetic  and  chemical  effects  in  real  systems.  Figure  4.40 
shows  the  calculated  solid  solubility  of  P  in  selected  alpha 
Fe-O.OIM-P  systems.  At  temperatures  above  the  system  Curie 
temperature,  the  solubility  curves  are  very  similar  to  those 
in  Figure  4.39.  However,  below  the  Curie  temperature,  the 
predicted  curves  deviate  greatly. 

Of  special  interest  is  the  Fe-V-P  system.  The  magnetic 
parameter  for  small  amounts  of  vanadium  should  be 
approximately  +1000  K.  A  positive  magnetic  parameter  raises 
the  solubility  limit  of  P.  Therefore,  although  the  Cr  curve 
depends  on  a  smaller  I'  ,  it  also  depends  on  a  smaller 
magnetic  parameter,  and  at  some  low  temperature  actually 
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crosses  below  the  V  curve.  Two  specific  examples  of  alpha 
Fe-M-P  systems  are  now  shown.  Dashed  lines  in  the  next  two 
figures  represent  the  paramagnetic  limits  of  solid 
solubility.  Figure  4.41  exhibits  the  predicted  alpha  Fe-Mo-P 
systems  and  Figure  4.42  the  predicted  Fe-V-P  system.  It  is 
clear  that  ferromagnetism  can  either  raise  or  lower  the 
solubility  depending  on  the  magnetic  parameter. 


4.12  Vacancy  Concentration  in  alpha  Fe-X 

The  following  treatment  for  ternary  systems  corresponds 
to  that  found  in  section  4.5  for  binaries.  That  is,  consider 
a  hypothetical  alpha  Fe-X  solid  solution,  denoted  A,  that 
contains  no  vacancies  at  any  temperature.  Now,  assume  that 
vacancies  introduced  into  the  solution  may  be  treated  like  a 
solute  species.  Therefore,  one  has  an  alpha  Fe-X-v  ternary 
solid  solution.  By  analogy  with  an  alpha  Fe-M-X  regular 
solid  solution,  the  free  energy  of  alpha  Fe-X-v  may  be 
written  as 

gT  =  XA°Gnm  +  Xv°Gnm  +  XX  Gnm  +  XAXX!nm  +  XXXvInm 

+  X.X  IAv  +  RT[X.£nX,  +  Xy£nXx  +  X  inX  ] 

Avnm  AA  aa  vv 

+  [l-HvXv][l-MxXx][G^g(TT)-[ATxXx+ATvXv]US^g]  (4143) 

The  equilibrium  vacancy  concentration  occurs  where  the 
ternary  free  energy  is  a  minimum  with  respect  to  the  vacancy 
concentration.  Again,  a  dilute  bulk  approximation  may  be 
applied,  XA  -  1 -X x,  and  the  free  energy  of  the  alpha  Fe-X 
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Figure  4.41 

The  predicted  solubility  of 
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Figure  4.42 

The  predicted  solubility  of 


solution  approximates  that  of  the  hypothetical  solution. 
Therefore,  the  derivative  may  be  written  as 


153 


T 


v 


=  0 


(4.144) 


where  the  magnetic  component  of  Q  is 


The  equilibrium  concentration  of  vacancies  is  given  by 


X  =  [1  -XM]exp(-[Q  +Q  ]/RT ) 
v  L  MJ  K  v  LVnm  ymaq J/  ' 


(4.146) 


nm  ^mag 


This  equation  is  similar  to  equation  (4.68).  A  large 
negative  magnetic  parameter  for  the  solute  X  slightly 
increases  the  vacancy  concentration.13 

If  the  solute  is  substitutional,  the  diffusion 
coefficient  of  X  in  alpha  iron  will  depend  on  the  vacancy 
concentration  in  the  bulk.  Therefore,  the  diffusion 
coefficient  of  X  will  behave  in  a  manner  similar  to  the 
sel f-di f fusion  coefficient  of  alpha  iron.  This  has  been 
observed  experimentally  for  various  species  in  Fe  [35,66]. 
The  definition  of  T  ^  is 


T 


T 


13There  is  a  solute-vacancy  attraction. 
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4.13  Grain  Boundary  Segregation  in  Alpha  Fe-M-X  Solid 
Solut ions 

The  thermodynamic  analysis  of  section  4.6  is  now 
applied  to  ternary  solutions  of  alpha  iron.  The  segregation 
free  energy  as  defined  in  equation  (4.85)  depends  on  the 
standard  free  energies  of  the  solute  components  and  their 
interaction  parameters.  Guttmann  [17]  proposed  that  the 
segregation  free  energies  for  an  alpha  Fe-M-X  system  may  be 
approximated  by 


AGX  =  °AGX  +  I  '  [X*-Xjj] 

AGM  =  °AGM  +  I'[X*-X“]  (4.14 
in  a  competitive  model.14  I7  is  a  preferential  interaction 
parameter  for  solute  X  and  M  such  that 


T  ,  _  tMX  TFeX  T  FeM 

1  "  1  '  1  1  (4.148) 

The  effect  of  an  attractive  I7  interaction  (positive  I7)  is 
to  increase  the  segregation  free  energy  and  therefore  the 
tendency  of  a  solute  to  segregate.  The  effect  of  I7  on 
segregation  free  energy  (in  absence  of  magnetic 
interactions)  is  shown  in  Figure  4.43  for  alpha  Fe-M-X.  The 
corresponding  calculated  grain  boundary  solute  concentration 
curves  for  competitive  segregation  are  given  in  Figure  4.44 
These  curves  were  produced  by  the  program  developed  in  the 
present  work  and  they  correspond  exactly  to  Guttmann7 s 


14  In  a  non-competitive  model,  I7  is  modified  by  the 
composition  of  an  M-X  compound. 
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Figure  4.44 

The  effect  of  I'  on  the  boundary  solute  content  f 
(A  =  20  KJ/mole,  B  =  40  KJ/mole,  C  =  60  K J / mo  1 e ) . 
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curves . 

The  assumption  that  the  segregation  free  energies  vary 
linearly  with  boundary  solute  content  should  allow  I#  and 
standard  segregation  free  energies  to  be  determined 
experimentally.  This,  in  fact,  has  been  investigated  by 
Guttmann  [60]  who  determined  that  the  standard  segregation 
free  energy  for  phosphorus  is  between  40  and  45  KJ/mole  and 
for  transition  metals  between  0  and  10  KJ/mole.  These  values 
will  be  used  here  for  convenience.  However,  as  is  shown  in 
the  present  analysis,  segregation  free  energy  does  not  vary 
linearly  with  boundary  solute  content.  This  is  evident  upon 
inspecting  the  magnetic  component  of  segregation  free 
energy.  This  affects  the  experimental  determinations  of  both 
I'  and  the  standard  free  energies.  Guttmann' s  free  energy  of 
segregation  is  approximately  equivalent  to  a  paramagnetic 
component  of  segregation  free  energy  and  is  used  as  such  in 
this  work. 

Guttmann' s  analysis  is  now  applied  to  ternary  alpha 
Fe-M-X  solid  solutions  using  the  chemical  potentials 
obtained  in  section  4.10.  Chemical  interactions  are  taken 
into  consideration  when  examples  of  specific  systems  are 
exhi bi ted . 

4.13.1  Case  I:  Grain  Boundary  Free  Energy  Without  Magnetic 
Terms 

This  approximation  is  considered  for  the  same  reasons 
stated  in  section  4.7.1,  i.e.  to  obtain  trends.  The  grain 
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boundaries  will  be  considered  non-magnet ic ,  i.e.  they  cannot 
take  part  in  any  magnetic  interactions  and  there  are  no 
magnetic  terms  in  the  boundary  free  energies.  Direct 
analytic  solutions  were  therefore  calculated.  To  obtain  an 
equilibrium  segregation  equation,  the  following  conditions 
are  applied  for  some  system  temperature  T, 


d>  X  aX 

=  y 

d>M  aM 

yy  =  y 

d>  F  e  aFe 

yy  =  y 


(4.149) 


The  chemical  potentials  from  section  4.10  are  defined, 
as  in  section  4.7.1,  to  be 


4>X  =  o  4>X  +  RT£)lXJ  .  ou 
M  Mnm  X  X 

■  °^nm  +  RUnXX  +  ^mag  <4.150) 

aX 

=  y 


y 


4>M 


o  d>M 
y 

Mnm 


RT£nXjjj 


o  aM 


+ 


RU*XM  + 


i  aM 
ymag 


aM 

y 


(4.151) 


d>Fe 

y 


o  d>Fe 

y 

pnm 

o  aFe 
^nm 


+  RT£nX 


+  RT£»iX 


<P 

Fe 

a 

Fe 


-  aw 


Fe 

aFe 

lmag 


aFe 

y 


(4.152) 
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Again,  by  eliminating  the  surface  tension  terms  and 
assuming  the  atomic  structure  of  the  grain  boundary  is  not 
altered  (all  partial  molar  areas  are  equal),  segregation 
equations  are  obtained, 


RT£n(X*X“e/X“X*e) 

RT£n(X*X“e/X“X^e) 


M 

mag 


o  aM  o  d>M  o  aFe^o  4>Fe,  aM 
u  -  u  -  u  +  u  +u 
Mnm  Mnm  Mnm  Mnm  Hmag 


4>Fe 

ymag 


(4.154) 


The  magnetic  segregation  free  energy  of  solute  X  is  given  by 


AG 


X 

mag 


-Mx[l-MMX“][G^g(T“)-[ATxx“+ATMX“]Us^] 

-  [l-MMX“][l-MxX“]ATx[S^g(T“)-V®g] 


The  segregation  free  energy  of  M  is  perfectly  symmetrical  to 
this  expression  (subscript  X  is  replaced  by  subscript  M,  and 
vice  versa).  If  the  amount  of  non-magnetic  solute  X  in  the 
bulk  is  small,  and  applying  the  approximation  in  equation 
(4.35),  then  the  segregation  free  energy  becomes 

AGmag  =  '  ^ 1 {Gmag  < T  > 

+  tATMXM‘ATX^Sinag(Ta)'  Smag^}  (4.156) 
Therefore,  unless  the  bulk  content  of  M  is  very  high,  the 
magnetic  segregation  free  energy  of  X  is  not  very  dependent 
on  the  magnetic  parameter  of  M.  An  example  of  this  situation 
is  shown  in  Figure  4.45.  These  curves  have  the  same  shapes 
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as  observed  in  the  analysis  of  binary  solid  solutions.  The 
competitive  boundary  segregation  curves  corresponding  to 
Figure  4.45  are  shown  in  Figure  4.46.  For  the  sake  of 
comparison,  the  bulk  solute  content  of  M  is  set  higher  than 
that  of  X,  simulating  the  situation  where  M  is  a  transition 
metal  alloying  element  and  X  is  an  impurity  element.  Due  to 
competition,  an  M  species  with  a  large  negative  magnetic 
parameter  will  depress  the  grain  boundary  content  of  X.  As 
the  magnetic  parameter  of  M  becomes  larger,  the  segregation 
curve  for  X  approaches  a  limit  which  corresponds  to  the 
equivalent  segregation  curve  in  the  alpha  Fe-X  binary. 

The  bulk  content  of  one  solute  may  also  drastically 
affect  the  grain  boundary  concentration  of  the  other. 

Figure  4.47  shows  the  behavior  of  predicted  grain  boundary 
solute  concentrations  with  varying  bulk  content  of  M.  A 
large  amount  of  solute  M  in  the  bulk  may  reduce  the  grain 
boundary  concentration  of  X  especially  at  higher 
temperature.  This  may  occur  even  if  M  has  little  effect  on 
the  magnetic  properties  of  the  solution  (magnetic  parameter 
of  M  is  zero).  This  situation  could  be  encountered  where  M 
is  a  metal  alloying  element  and  X  an  impurity  element. 

Case  I,  as  outlined,  points  out  the  general  trends  of 
the  effect  of  ferromagnetism  on  segregation.  However,  the 
shortcomings  of  Case  I  discussed  for  binary  solutions  are 
likewise  present  here.  Therefore,  Case  II  will  now  be 
considered . 
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Figure  4.46 

The  effect  of  the  magnetic  parameter  of  M  on  boundary  solute  contents  (A  =  -1000 
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4.13.2  Case  II:  Grain  Boundaries  May  Be  Ferromagnet ic 

As  described  in  section  4.7.2,  the  only  difference 
between  a  2-d  and  a  3-d  phase  is  the  2-d  surface  tension. 

The  chemical  potentials  of  the  grain  boundary  species  would 
have  magnetic  components  and  would  be  written  analogously  to 
the  chemical  potentials  in  section  4.10.  The  Guttmann 
analysis  is  again  applied15  and  for  a  system  temperature  T, 
the  equilibrium  segregation  equations  result,  from  the 
cond i t i ons 


y 


<t>i 


o  <j)i 


y  +  R  UhXT 
Mnm  l 


. 


CFO  . 


y 


<f>i 

mag 


°ual  +  RUkX?  +  y 


nm 


ai 

mag 


ai 

y 


.(4. 157) 


The  magnetic  segregation  of  solute  X  is  given  by 

RU»(X*X“e/X“X*e)  -  AG*m  +  nAG^g 

o  aX  o  <|>X  o  aFe  .  o  4>Fe 

=  u  u  -  u  +  u 

Hnm  nm  nm  nm 

aX  aFe  X  4>Fe 

+  u  -y  -  y  +y 

Mmag  Mmag  Mmag  Mmag  (4.158) 

The  magnetic  free  energy  of  segregation  for  solute  M  is 
symmetrical  to  the  above  equation  with  subscript  X  being 
replaced  by  M  and  vice  versa.  After  treatment  as  in  Case  I, 


15The  partial  molar  areas  for  all  species  are  equal. 
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n.j 

AG 

mag 


MX[1-V^t<eag(T^-^TXXX  +  ATMX^Usmag] 

'  HXCl-MHXS][Gmag(Ta)-[ATXXX+iTMXM]Usmag^ 

-  [l-MMX*][l-HxX$]ATx[S^ag(T^)-uS^g] 

-  n-MMX“][l-HxX“]ATx[S^g(Ta)-US^g] 


(4.159) 


This  equation  cannot  be  evaluated  directly  but  equation 


(4.85)  must  be  considered  for  X  and  M,  such  that 


xj  =  X“  exp ( 1 1 AGX/RT ) 

1'XX'XM+XXexp( ilAGX/RT)+XMexp( IIagM/RT) 

X$  =  X“  exp(HAGM/RT) 

1  'Xx’XM+XXexP(  UAGX/RT)  +  X“exp(  nAGM/RT) 

x  n  x  M  (4.160) 


are  solved  simultaneously  for  roots.  This  is  a  somewhat  more 
difficult  numerical  analysis  than  performed  in  the  binary 
case . 


Figure  4.48  shows  the  magnetic  segregation  free 
energies  of  alpha  Fe-M-X  as  a  function  of  temperature  and 
magnetic  parameter  of  M.  As  in  Case  I,  large  negative 
magnetic  parameters  lead  to  larger  segregation  free 
energies.  The  general  shapes  of  the  binary  curves  are 
retained.  However,  since  X  and  M  interact  and  are 
competitive,  a  reduced  magnetic  segregation  free  energy  for 
M  (due  to  a  larger  magnetic  parameter  of  M)  reduces  the 
magnetic  segregation  free  energy  for  X.  These  curves  may 
have  steps  and  a  complex  behavior  (due  to  magnetic 


166 


”  0001 


00M 


0021 


008 


009 


A 

X 

KS 


UJ 

a. 

D 

h- 

< 

£ 

UJ 

Q. 

Z 

UJ 

H 


00* 


002 


0 


□ 

G 

Q 

G 

G 

G 

G 

G 

g 

G 

G 

G 

G 

G 

G 

□ 

Q 

G 

G 

G 

G 

G 

□ 

ID 

G 

in 

G 

in 

in 

(T) 

CM 

CM 

H 

H 

l 

C3“10W/r>  *03S  30  A0d3N3  33«d  *OVW 


Figure  4.48 

The  effect  of  the  magnetic  parameter  of  M  on  the  magnetic  segregation  free  enerqy 
according  to  Case  II  (A  =  -1000  K,  B  =  0  K,  C  =  1000  K). 
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interaction  and  site  competition).  Figure  4.49  is  the 
corresponding  grain  boundary  solute  content  curve.  Curves  A 
show  that  similar  solute  share  grain  boundary  sites.  A  large 
difference  in  segregation  free  energies  usually  allows  one 
solute  species  to  'win  out'  over  the  other.  The  effect  of 
bulk  content  of  M  on  the  magnetic  segregation  free  energies 
is  shown  in  Figure  4.50  and  the  corresponding  grain  boundary 
segregation  curves  in  Figure  4.51.  It  is  evident  that  in  a 
competitive  model  even  solute  M  with  a  zero  magnetic 
parameter  but  large  bulk  content  may  reduce  the  boundary 
content  of  X.  A  similar  effect  from  the  non-magnetic 
segregation  free  energy  of  M  is  shown  in  Figure  4.52. 

In  a  noncompetitive  situation,  the  grain  boundary 
segregation  curves  are  similar  to  the  binary  segregation 
curves.  That  is,  the  segregation  curve  for  each  solute  is 
shifted  when  the  segregation  free  energy  for  that  solute  is 
modified  by  considering  magnetic  terms.  The  effect  of  one 
solute  on  the  other  is  small  unless  there  is  a  strong 
(usually  chemical)  interaction  between  the  solute.  In  such  a 
case,  the  strongly  segregating  solute  'drags'  the  other 
solute  to  the  boundary,  so  that  the  segregation  of  the 
weakly  segregating  solute  is  enhanced. 
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Figure  4.49 

Predicted  grain  boundary  solute  contents  for  alpha  Fe-O.OIM-X  according  to  Case  II 
(A  =  -1000  K,  B  =  0  K,  C  =  1000  K) . 
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Figure  4.50 

The  effect  of  bulk  content  of  solute  M  on  the  magnetic  segregation  free  enerqies 
according  to  Case  II  (A  =  0.1,  B  =  0.01,  C  =  0.001). 
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Figure  4.52 

The  effect  of  the  non-magnet ic  segregation  free  energy  of  M  on  grain  boundary 
solute  content  acccording  to  Case  II  (A  =  50  KJ/mole,  B  =  40  KJ/mole, 

C  =  30  KJ/mole ) . 
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4.14  Predicted  Grain  Boundary  States  in  Certain  Ternary 
Systems 

As  outlined  in  section  4.8,  the  state  of  a  grain 
boundary  is  dependent  on  its  solute  content.  The  same 
possibility  exists  that  boundaries  rich  in  solute  can 
achieve  a  transition  to  the  paramagnetic  state  as 
segregation  causes  the  Curie  temperature  of  the  boundary  to 
fall  below  the  system  temperature.  The  grain  boundary  Curie 
temperature  for  an  alpha  Fe-M-X  is 


T 


<p  _  T  F  e 


<P 

XAX 


+  AT  v  X  „  +  AT..X 


<P 
M'*M 


(4.161) 


No  simple  criteria  for  transition  may  be  applied  to  ternary 
segregation  curves  as  was  the  case  in  the  binary  model.  This 
is  due  to  the  complex  relationship  between  the  various  grain 
boundary  solute  concentrations. 

It  is  possible,  however,  to  represent  the  transitions 
as  was  done  for  binaries  in  Figure  4.33.  Figure  4.53  shows 
the  equilibrium  grain  boundary  segregation  curves  equivalent 
to  the  segregation  curves  in  Figure  4.52.  With  lowering 
temperatures,  the  paramagnetic  to  ferromagnetic  to 
paramagnetic  transitions  are  easily  seen. 

Three  specific  competitive  alpha  Fe-.OIM-P  systems  are 
now  considered.  Figure  4.54  shows  the  predicted  grain 
boundary  concentration  of  P  for  alpha  Fe- . 0 1 Ni - . 000 1 P  under 
chemical,  magnetic  or  both  types  of  interactions.  As  is 
irmedi ately  evident,  the  magnetic  interaction  significantly 
shifts  the  segregation  curve  to  higher  temperatures 
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(increases  segregation).  In  this  situation,  the  magnetic 
interaction  is  stronger  than  the  chemical  interaction. 

The  second  example  concerns  alpha  Fe- . 0 1 Mn- . 000 1 P . 
Figure  4.55  shows  the  predicted  grain  boundary  concentration 
of  P  under  the  various  interaction  possibilities.  The 
magnetic  and  chemical  interactions  are  of  the  same  order  of 
strength.  The  combination  of  both  interactions  shifts  the 
segregation  to  significantly  higher  temperatures. 

However,  the  third  example  shows  an  alpha 
Fe- . 0 1 Cr - . 000 1 P  system.  Figure  4.56  indicates  that  the 
excess  chemical  segregation  free  energy  has  more  impact  on 
the  grain  boundary  composition  than  the  magnetic  segregation 
free  energy.  The  calculated  segregation  curve  is  still 
significantly  shifted  to  higher  temperatures  (dotted  line) 
but  the  step  or  transition  temperature  is  not  as  affected. 
Figure  4.57  shows  the  corresponding  equilibrium  Curie 
temperature  curves.  This  indicates  that  strong  chemical 
effects  may  also  reduce  the  range  over  which  the  grain 
boundary  may  be  ferromagnetic. 

A  word  must  be  said  about  the  steps  occurring  in  the 
previous  figure.  When  the  segregation  curve  shows  a 
non-function  relation,  that  is,  when  there  is  more  than  one 
grain  boundary  composition  possible  for  a  system 
temperature,  the  solid  solution  becomes  metastable.  The 
stable  configuration  requires  two  solid  solutions  with  a 
transition  occurring  at  the  inflection  in  the  segregation 
curve.  This  point  occurs  at  the  mid-point  of  the  S-shaped 
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Figure  4.57 

Equilibrium  Curie  temperature  curves  for  alpha  Fe-Cr-P  corresponding  to  Fiqure  4.56 
(C  -  chemical  only,  M  =  magnetic  only,  B  =  chemical  and  magnetic). 
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curve.  However,  for  the  magnetic  segregation  free  energy, 
the  inflection  and  transition  do  not  occur  at  the  mid-point 
but  rather  depend  on  equation  (4.161).  These  steps  in  the 
previous  figures  represent  transtions  from  one  state  to 
another.  When  both  chemical  and  magnetic  interactions  are 
considered,  the  temperature  of  the  inflection  point  seems 
dominated  by  equation  (4.161).  Therefore,  even  if  the 
calculated  segregation  curve  shifts  a  great  deal  to  the 
right,  the  transition  temperature  may  not  shift 
significantly  as  is  the  case  in  the  Fe-Cr-P  system  (Figure 
4.56)  . 

4.15  Non-linear  Changes  in  Curie  Temperature 

The  assumption  that  the  Curie  temperature  of  iron 
varies  linearly  with  solute  addition  is,  of  course,  only  a 
first  order  approximation.  A  simple  glance  at  any 
compilation  of  binary  phase  diagrams  shows  examples  of 
non-linear  variations  of  Curie  temperature.  In  iron,  this 
property  is  especially  noticeable  where  a  solute  has  a  large 
solid  solubility.  Examples  include  Fe-Cr,  Fe-Co  and  Fe-V. 

The  Fe-V  is  especially  interesting  since  the  Curie 
temperature  rises  rapidly  to  about  10  at%  V  and  then  drops. 

Figure  4.58  shows  the  grain  boundary  content  of  a 
solute  X  wi th  a  linear  and  non-linear  magnetic  parameter. 

The  linear  magnetic  parameter  is  1000  K  while  the  non-linear 
magnetic  parameter  is  approximated  by  using  1000  K  until 
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10  at%  X  and  -1270  K  thereafter.  This  example  gives  an  idea 
of  the  complex  segregation  behavior  possible  with  non-linear 
magnetic  parameters,  specifically  V. 


4.16  A  Calculated  Magnetic  Parameter 

It  is  possible  to  theoretically  calculate  the  magnetic 
parameters  for  iron  using  a  statistical  mechanical  model  of 
fer romagnet i c  binary  alloys  containing  a  non-magnetic 
solute.  Arita  [67]  has  presented  a  quasi -chemica 1  model  of 
ferromagnetic  binary  alloys  using  a  strictly  regular  solid 
solution.  His  model  is  presented  here  because  of  its 
relevance  to  the  present  work.  The  thermodynamic  properties 
of  a  ferromagnet ic  alloy  were  described  by  the  same  methods 
as  may  be  applied  to  non-magnetic  mixtures.  The 
discontinuity  in  the  heat  capacity  of  the  alloy  was 
attributed  to  the  Weiss  field  (exchange  interaction  between 
the  spins  of  the  electrons  of  neighbouring  atoms)  and  the 
exchange  interaction  was  compared  to  the  change  of  the 
pair-wise  interaction  energy  upon  mixing  of  two  different 
atomic  species.  The  basic  assumptions  of  Arita' s  model  are 
the  same  as  presented  in  this  work,  that  is, 

1.  there  is  only  one  type  of  site  available  for  all  atoms 
present , 

2.  the  atomic  volumes  of  the  species  are  alike  and 
unaltered  on  mixing, 

3.  only  nearest  neighbour  interactions  are  considered,  and 
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4.  the  vibrational  partition  function  of  an  atom  is 
unchanged  upon  mixing. 

Briefly,  a  partition  function  (P.F.)  was  defined  on 
strictly  regular  behavior  and  the  above  assumptions.  A  free 
energy  function  for  the  binary  solid  solution  was  obtained 
by 

G  -  F  =  -  kT-MP.F.) 

(4.162 

A  free  energy  of  mixing  could  then  be  calculated.  The  free 
energy  of  the  binary  was  mathematically  manipulated  to 
obtain  the  initial  slope  of  the  Curie  temperature  with 
addition  of  a  non-magnetic  solute.  Unfortunately,  the  Curie 
temperature  function  could  not  be  analytically  determined, 
however  the  slope  was  given  as 


(dT  /dX  ' 
v  c 


=  -4kT2  exp(-4w/kT  ) 

C 


(4. 163) 


where 


j  is  the  positive  exchange  integral, 

z  is  the  number  of  nearest  neighbours, 

w  ..  is  the  interaction  energy  of  a  pair  of  i  and  h 


ih 


atoms,  and 


w  -  w1 2  -  [wn  -w22] 


2 


(4. 164) 


which  is  essentially  the  same  as  the  magnetic  parameter  used 
in  the  present  work.  The  parameter  w  may  be  determined  from 
the  free  energy  of  mixing  and  experimental  thermodynamic 
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data  such  as  [68].  Little  data  is  available  on  the  initial 
slope  of  the  Curie  temperature  for  iron  alloys;  however,  the 
4  data  points  noted  by  Arita  would  fit  in  well  with  the 
magnetic  parameters  used  in  this  work. 


4.17  Applicability  of  the  model 

The  equations  used  to  model  the  equilibrium  segregation 
behavior  in  fer romagnet i c  iron  actually  describe  this 
behavior  for  similar  iron  grains  separated  by  planar 
boundaries  containing  only  one  type  of  site  for  the 
competitive  situation  and  two  types  of  sites  for  the 
non-competitive  ternary  situation.  Although  all  grains  have 
a  definite  solute  concentration  there  were  no  bounds  placed 
on  the  system. 

There  was  no  information  given  about  the  'size'  of  the 
iron  grains  (and  therefore  about  the  'size'  of  the 
boundary).  The  model  reacts  as  if  the  iron  crystals  are 
infinitely  large  and  there  is  always  enough  solute  present 
to  saturate  the  boundary.  This  situation  occurs  when  the 
exponential  terms  in  the  segregation  equations  dominate  the 
bulk  solute  content  (as  happens  at  very  low  temperatures ) . 
Therefore,  to  include  solid  solubility  effects  at  lower 
temperatures  in  the  segregation  model,  some  dimensional 
information  must  be  provided.  In  general,  the  present  model 
should  be  valid  above  room  temperature. 
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All  thermodynamic  functions  used  were  obtained  from 
experimental  and  theoretical  data  for  polycrystalline 
iron16.  In  the  present  description  one  solute  species  may 
completely  displace  all  other  species  from  the  boundary  in  a 
competitive  situation.  This  occurrence  is  clearly 
unreasonable.  The  distribution  of  sites  in  the  boundary  is 
orientation  and  structure  dependent.  Although  any  type  of 
atomic  species  present  may  attempt  to  fill  a  particular 
site,  the  probability  of  a  particular  species  remaining  in  a 
particular  boundary  site  depends  on  the  relative  decrease 
(or  increase)  of  the  system  free  energy  by  the  binding 
energy  for  that  atom  in  that  site.  Competition  between 
solute  species  would  depend  on  their  relative  binding 
energies  for  the  sites  available.  Similar  binding  energies 
indicate  a  highly  competitive  situation  while  highly 
dissimilar  binding  energies  indicate  a  non-competitive 
si tuat ion . 

No  attempt  has  been  made  to  apply  the  present  work  to 
segregation  data  in  the  literature.  There  are  a  number  of 
reasons  for  this.  The  experimental  data  (see  chapter  2)  is 
generally  scattered  with  greatly  varying  saturation  levels 
and  generally  no  information  about  the  grain  boundaries 
analysed.  Also,  the  interpretation  of  AES  data  has  changed 
over  the  years  making  it  difficult  to  compare  data  from 
various  sources. 


1 6These  functions  are  thus  statistically  average  functions. 
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Finally,  various  workers  [69]  are  seeking  evidence  of 
two-dimensional  phase  transitions  at  grain  boundaries. 
Unfortunately,  there  have  been  no  conclusive  results 
publ i shed . 


5.  CONCLUSIONS 


The  conclusions  derived  during  the  course  of  this  work 

are  listed  below  in  point  form: 

1.  The  experimental  techniques  developed  were  successful  in 
producing  'pure'  iron  phosphides  and  alpha  Fe-P  alloys. 

2.  Using  EPMA  to  determine  the  solute  content  from  a 
precipitate  surface  into  the  bulk  provided  an  effective 
means  of  estimating  the  equilibrium  limit  of  solid 
solubility  of  the  solute. 

3.  Below  the  Curie  temperature,  an  anomaly  in  the  solid 
solubility  of  P  in  alpha  Fe  was  experimentally 

conf i rmed . 

4.  The  theory  of  Nishizawa,  applied  to  the  solid  solubility 
of  P  in  alpha  Fe,  resulted  in  a  calculated  magnetic 
parameter  of  -1200  K  for  P  confirming  Nishizawa' s 
magnetic  parameter. 

5.  An  anomaly  in  the  vacancy  concentration  of  iron  below 
the  Curie  temperature  was  predicted.  This  prediction  was 
consistent  with  sel f-di f fusion  data  for  iron. 

6.  The  ferromagnetic  behavior  of  iron  modifies  the  grain 
boundary  solute  segregation  behavior  for  binary  and 
ternary  iron-based  alloys  when  compared  to  purely 
non-magnet ic  behavior. 

7.  For  iron-based  binaries  a  non-magnetic  solute  which 
drastically  lowers  the  system  Curie  temperature  will 
also  'saturate'  grain  boundary  sites  at  a  higher 
temperature.  That  is,  the  segregation  curve  is  shifted 
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to  the  right  when  compared  to  an  analysis  which  neglects 
fer romagnet i sm. 

8.  Decreasing  solute  magnetic  parameters  increase  solute 
segregation  while  increasing  magnetic  parameters 
decrease  solute  segregation. 

9.  The  analytic  Case  I  approximations  successfully  show  the 
above  trends  for  alpha  Fe-X  binaries. 

10.  The  numerical  solutions  of  Case  II  for  alpha  Fe-X  also 
successfully  show  the  above  trends  but  do  so  in  a  more 
realistic  manner  with  no  magnetic  contributions  above 
the  system  Curie  temperature. 

11.  Two-dimensional  phase  transitions  are  predicted  to  occur 
in  the  grain  boundaries  of  alpha  Fe-X  binaries.  These 
transitions  occur  at  the  grain  boundary  Curie 
temperature  suggesting  magnetic  transitions. 

12.  Below  the  system  Curie  temperature  the  solubility  of  P 
in  alpha  Fe-M-P  deviates  from  the  paramagnetic  behavior 
in  a  manner  dependent  on  the  magnetic  parameters  of  M 
and  P.  Small  (negative)  magnetic  parameters  decrease  P 
solubility  while  large  (positive)  magnetic  parameters 
increase  P  solubility  as  compared  to  paramagnetic 
behavior . 

13.  The  vacancy  concentration  in  alpha  Fe-X  was  calculated. 

14.  In  iron-based  ternaries,  a  large  negative  magnetic 
parameter  for  one  of  the  solute  species  decreases  the 
segregation  of  the  other  solute. 

15.  Likewise,  increasing  the  bulk  content  of  one  solute 
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decreases  the  segregation  of  the  other. 

16.  The  segregation  behavior  of  a  solute  in  an  iron-based 
ternary  depends  on  the  relative  magnitudes  of  the 
magnetic  and  non-magnet ic  segregation  free  energies  and 
the  relative  positions  of  the  two-dimensional 

trans i t ions . 

17.  A  positive  magnetic  free  energy  for  a  solute  increases 
the  grain  boundary  segregation  of  that  solute. 

18.  The  two-dimensional  transitions  for  magnetic  systems 
occur  when  the  system  temperature  equals  the  grain 
boundary  Curie  temperature. 

19.  When  chemical  interaction  parameters  and  magnetic 
parameters  for  solute  are  Known,  the  segregation 
behavior  of  alpha  Fe-M-X  ternaries  may  be  predicted. 

20.  Non-linear  changes  in  the  system  Curie  temperature 
produce  complex  segregation  behavior. 


—  N 


6.  FUTURE  WORK 

There  are  a  myriad  of  possibilities  for  expanding  on 

the  present  research.  A  few  proposals  are  discussed  below: 

1.  Most  AES  segregation  data  is  obtained  from  complex 
polycrystalline  alloys  that  will  fracture 

i ntergranu 1 ar ly .  The  boundary  types  and  compositions  may 
vary  widely.  It  is  suggested  that  Fe-X  bicrystals  (or 
very  simple  polycrystals)  be  used  to  study  grain 
boundary  segregation  in  an  effort  to  control  the 
boundary  type  and  hopefully  the  boundary  sites 
available.  This  would  be  a  difficult  problem  to  solve  as 
a  brittle  fracture  must  be  made  to  occur  at  the  grain 
boundary  for  AES  analysis. 

2.  The  effect  of  any  other  2nd-order  transition  phenomena, 
e.g.  order -di sorder  phenomena,  on  the  grain  boundary 
solute  segregation  may  be  thermodynamically  modeled  in  a 
simi 1 ar  mannner . 

3.  The  present  analysis  may  be  modified  to  include  a 
spectrum  of  binding  energies  approach  when  more  data 
becomes  ava i 1 ab 1 e . 

4.  The  magnetic  properties  of  an  iron  single  crystal  are 
anisotropic.  There  may  be  an  effect  on  grain  boundary 
segregation  or  precipitate  morphology.  In  particular, 
precipitate  morphology  may  be  different  above  and  below 
the  system  Curie  temperature. 

5.  As  grain  boundaries  'saturated'  with  non-magnetic 
impurities  should  disrupt  magnetic  order  at  the 
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boundary,  it  should  be  interesting  to  observe  domain 
boundaries  that  cross  grain  boundaries  in  segregated  and 
unsegregated  iron  alloys. 

6.  Any  external  conditions,  e.g.  strong  magnetic  fields, 
that  affect  the  magnetic  behavior  of  an  iron  alloy  may 
also  affect  the  segregation  of  that  alloy. 

7.  The  predicted  anomaly  of  P  solubility  in  alpha  Fe-M-P 
should  be  experimentally  verified,  with  particular 
interest  in  the  Fe-V-P  system. 

8.  A  very  interesting  possibility  concerns  the  detection  of 
the  state  of  the  grain  boundary  ,e.g.  2-d  phase 
transitions.  A  great  deal  of  information  may  be  obtained 
from  the  spins  of  the  electrons  at  an  interface.  This 
would  require  the  use  of  a  polarized  electron  beam  in 
electron-optical  devices.  At  present  the  author  Knows  of 
no  such  device. 


X 

. 


191 


BIBLIOGRAPHY 


1.  C.L.  Briant  and  S.K.  B a ner j i , In tergranu 1 ar  Failure  in 
Steel  -  The  Role  of  Gra i n  Boundary  Compos i t  i ons ,  Repo r t 
No.  77CRD118,  General  Electric  Company  (1977). 

2.  J.H.  Westbrook  and  K.  Aust,  Acta  Met.  11,  1151  (1963). 

3.  H.  Wiedersich  and  S.  R.  MacEwen,  d.  Nucl .  Mater.  90,  157 
( 1980) . 

4.  E.  D.  Hondros  and  M.  P.  Seah,  Int.  Metal  1 .  Rev.  22,  262 
( 1977)  . 

5.  M.  Guttmann,  Met.  Trans.  8A,  1383  (1977). 

6.  R.  M.  Bozorth,  Ferromagnetism ,  Van  Nostrand  Co.,  New 
York  ( 1959) . 

7.  S.  Chikazumi ,  Physics  of  Magnetism,  Wiley,  New  York 
( 1964) . 

8.  A.  H.  Morrish,  The  Physical  Principles  of  Magnetism, 
Wiley,  New  York  (1965). 

9.  M.  C.  Inman  and  H.  R.  Tipler,  Metal  1 .  Rev.  8,  105 
( 1963)  . 

10.  J.  W.  Gibbs,  Collected  Works,  Vol.  1,  p.  219,  Yale  Univ. 
Press,  New  Haven  (1957). 

11.  A.  H.  Cottrell,  in  Report  of  a  conferenceon  the  strength 
of  solids,  p.  30,  Physical  Society,  London  (1948). 

12.  D.  McLean,  Grain  Boundaries  in  Metals,  Clarendon  Press, 
Oxford  (  1 957  )  . 

13.  M.  P.  Seah  and  E.  D.  Hondros,  Proc.  R.  Soc.  335A,  191 
( 1973)  . 

14.  S.  Baynard,  Master ' s  Thesi s ,  Ecole  Polytechnique  (1977). 

15.  S.  Brunauer  et  al.,  d.  Amer .  Chem.  Soc.  60,  309  (1938). 

16.  R.  H.  Fowler  and  E.  A.  Guggenheim,  Statistical 
Thermodynamics,  University  Press,  Cambridge  (1960). 

17.  M.  Guttmann,  Surf.  Sci .  53,  213  (1975). 


, 


192 


18.  M.  C.  Inman  and  H.  R.  Tipler,  Acta  Met.  6,  73  (1958). 

19.  P.  V.  Ramasubramani an  and  D.  F.  Stein,  Met.  Trans.  4, 
1735  (1973). 

20.  M.  Hansen,  Const itut ion  of  Binary  Alloys,  2  ed., 
McGraw-Hill,  New  York  (1958).  -  R.  P.  Elliott,  First 
supplement  (1965).  -  F.  A.  Shunk,  Second  supplement 

( 1969)  . 

21.  A.  E.  Vol ,  Handbook  of  Binary  Metallic  Sysy terns; 
Structure  and  Properties  (translated  by  J.  Schmorak), 
Israel  Program  for  Scientific  Translation,  Jerusalem 

( 1966)  . 

22.  J.  L.  Haughton,  J.  Iron  Steel  Inst.  119,  417  (1927). 

23.  H.  Kaneko,  T.  Nishizawa,  K.  Tamaki  and  A.  Tanifuji,  d. 
Japan  Inst.  Metals  29,  166  (1965). 

24.  E.  Hornbogen,  Trans.  Amer .  Soc.  Metals  53,  569 

25.  E.  Hornbogen,  Trans.  Amer.  Soc.  Metals  55,  719 

26.  E.  Hornbogen,  Trans.  Amer.  Soc.  Metals  56,  16 

27.  H.  P.  Hoffmann,  K  Lohberg  and  W.  Reif,  Arch 

E i senhuttenw .  41,  975  (1970). 

28.  A.  S.  Doan,  Jr.  and  J.  I.  Goldstein,  Met.  Trans.  1,  1759 
( 1970)  . 

29.  M.  Ko  and  T.  Nishizawa,  J .  Japan  Inst.  Met.  43,  188 
( 1979)  . 

30.  W.  D.  Jones,  Fundamental  Principles  of  Powder 
Metallurgy,  E.  Arnold,  London  (1960). 

31.  P.  Lindskag,  J.  Tengzelius  and  S.  A.  Kvist,  Modern 
Developments  in  Powder  Metal  1 urgy  10,  97  (1976). 

32.  G.  Seibel ,  Mem.  Sci .  Rev.  Met.  61,  413  (1964). 

33.  K.  C.  Russell  and  H.  I.  Aaronson,  J.  Mater.  Sci.  10, 

1991  (1975). 

34.  J.  R.  Bradley  and  H.  I.  Aaronson,  Met.  Trans.  8A,  317 
( 1977)  . 

35.  G.  Luckman,  R.  A.  Didio  and  W.  R.  Graham,  Met.  Trans. 
12A,  253  (1981). 


193 


36.  S.  D.  Conte,  Elementary  Numerical  Analysis ,  chapter  2, 
McGraw-Hill,  Toronto  (1965). 

37.  R.  G.  E.  Franks,  Modeling  and  Simulation  in  Chemical 
Engineering,  p.  26,  Wi ley- Interscience ,  Toronto  (1972). 

38.  K.  E.  Szklarz  and  M.  L.  Wayman,  Acta  Met.  29,  341 
(1981). 

39.  C.  Zener,  d.  Metals  7,  619  (1955). 

40.  R.  J.  Weiss  and  K.  J.  Tauer,  Phys.  Rev.  103,  1490 
(  1956) . 

41.  R.  J.  Weiss  and  K.  J.  Tauer,  d.  Phys.  and  Chem.  Solids 
4,  135  (1958). 

42.  A.  P.  Miodownik,  Cal phad  1,  133  (1977). 

43.  G.  Inden,  Z.  Metal  IK.  66,  577  (1975). 

44.  T.  Nishizawa,  M.  Hasebe  and  M.  Ko,  Acta  Met.  27,  817 
(  1979)  . 

45.  M.  Hillert,  T.  Wada  and  H.  Wada,  d.  Iron  Steel  Inst. 

205,  539  (  1967)  . 

46.  H.  Harvig,  G.  Kirchnar  and  M.  Hillert,  Met.  Trans.  3, 

329  ( 1972) . 

47.  Metals  Handbook,  8th  ed . ,  Vol .  8,  American  Society  for 
Metal s . 

48.  T.  Takayama,  M.  Y.  Wey  and  T.  Nishizawa,  d.  dapan  Inst. 
Metals  45,  (1981)  in  press 

49.  G.  R.  Speich,  J.  A.  Gula  and  R.  M.  Fisher,  in  The 
Electron  Microprobe  (Edited  by  T.  D.  McKinley,  K.  F.  J. 
Heinrich  and  D.  B.  Witty),  p.  525,  Wiley,  New  York 

( 1966)  . 

50.  A.  H.  Qureshi ,  Z.  MetallK.  52,  799  (1961). 

51.  H.  A.  Wriedt  and  L.  S.  Darken,  Trans.  Met.  Soc.  A.I.M.E. 
218,  30  ( 1960) . 

52.  G.  R.  Speich,  L.  Zwell  and  H.  A.  Wriedt,  Trans.  Met. 

Soc.  A.I.M.E.  230,  939  (1964). 

53.  R.  J.  Borg  and  C.  E.  Birchenall,  Trans.  Met.  Soc. 

A.I.M.E  218,  980  (  1960) . 


' 

h&  1 1  i  i(  .*> 


194 


54.  F.  S.  Buffington,  K.  Hirano  and  M.  Cohen,  Acta  Met.  9, 
434  ( 1961 ) . 

55.  C.  M.  Walter  and  N.  L.  Peterson,  Phys.  Rev.  178,  922 
(1968). 

56.  G.  Hettich,  H.  Meher  and  K.  Maier,  Scripta  Met.  11,  795 
(  1977) . 

57.  H.  E.  Schaefer,  K.  Maier,  M.  Weller,  D.  Herlach,  A. 
Seeger  and  J.  Diehl,  Scripta  Met.  11,  803  (1977). 

58.  M.  Guttmann  and  D.  McLean,  in  Interfacial  Segregation, 
p.  261,  American  Society  for  Metals  (1979). 

59.  R.  Defay  and  I.  Prigogine,  Surface  Tension  and 
Adsorption  (Translated  by  D.  H.  Everett),  Wiley,  New 
York  (  1966)  . 

60.  M.  Guttmann,  private  communication 

61.  C.  L.  White  and  W.  A.  Coghlan,  Met.  Trans.  8A,  1403 
(  1977)  . 

62.  R.  J.  Harrison,  G.  A.  Bruggeman  and  G.  H.  Bishop,  in 
Grain  Boundaries:  Structure  and  Properties,  G.  A. 
Chadwick  and  D.  A.  Smith,  eds . ,  p.  45,  Academic  Press, 
New  York  (  1976) . 

63.  H.  G lei  ter  and  B.  Chalmers,  Prog.  Mat.  Sci .  16  (B. 
Chalmers,  J.  W.  Christian  and  T.  B.  Massalski,  eds.), 
Pergamon  Press,  Oxford  (1972). 

64.  M.  Guttmann,  Met.  Sci.  10,  337  (1976). 

65.  M.  Hillert,  in  Phase  Transformations,  p.  181,  American 
Society  for  Metals  (1970). 

66.  R.  J.  Borg  and  D.  Y.  F.  Lai,  Acta  Met .  11,  861  (1963). 

67.  M.  Arita,  Acta  Met .  26,  259  (1978). 

68.  R.  Hultgren  et  a  1 . ,  Selected  Values  of  the  Thermodynamic 

Properties  of  Binary  Alloys,  American  Society  for  Metals 
(  1973) . 

69.  W.  Losch  and  I.  Andreoni ,  Scripta  Met.  12,  277  (1978). 


. 


195 


APPENDIX  1 

This  appendix  contains  a  list  of  the  major  symbols  used  in 
this  work. 

subscript  'mag'  indicates  a  magnetic  component 
subscript  ' nm'  indicates  a  non-magnetic  component 
a  =  atomic  fraction  of  M  in  M  X 
a  =  act i vi ty  of  i 

A  =  1.08444  J/mole/K  in  equation  (4.9) 

b  =  atomic  fraction  of  X  in  M  X 

B  =  0.67432  J/mole/K  in  equation  (4.9) 

C  =  a  constant  in  Q  =  C  +  DT,  d/mole 

Cp1  =  heat  capacity  at  constant  pressure  of  i 

D  =  a  constant  in.  Q  =  C  +  DT,  J/mole/K 

Dj  =  diffusion  coeficient  of  i  in  alpha  iron 

E  =  segregation  binding  energy 

F  1  =  Helmholtz  free  energy  of  i 

G1  =  Gibbs  free  energy  of  i 

0  G 1  =  standard  free  energy  of  i 

H  1  =  enthalpy  of  i 
•  • 

I =  interaction  parameter  for  i  and  j 
r  =  Guttmann  interaction  parameter 
j  =  exchange  integral 
k  =  Boltzmann's  constant 

L  =  system  temperature  divided  by  Curie  temperature 
m  =  4.766  in  equation  (4.9) 

M  =  solute  M  or  a  constant  M  in  the  binary  model 
M -j  =  solute  intensity  factor  of  i 


N.j  =  number  of  atoms  of  i 
n  =  3.613  in  equation  (4.9) 

P  =  pressure  on  system 
Q  =  free  energy  of  solution  for  X  in  Fe 
Q  =  free  energy  of  vacancy  formation  in  Fe 
=  gas  constant 
S1  =  entropy  of  i 
T  =  system  temperature 
T1  =  adjusted  temperature  for  i 
T c  =  Curie  temperature  for  i 
v  =  vacancy 

V  =  system  volume 

=  exchange  energy  between  atoms  i  and 
X  -  solute  X 

X  =  atomic  fraction  of  i 
X  =  limit  of  solubility  of  X  in  Fe 

Y  =  atomic  fraction  of  sites  filled  by  i 
z  =  coordination  number 

a  =  bulk  phase 
4>  =  boundary  phase 

aG1  =  free  energy  of  segregation  for  i 
AGm  =  free  energy  of  mixing 
aT.|  =  magnetic  parameter  of  i 

e  =  ei  =  Feaxb 

62  =  MaXb 
■ 

y1  =  chemical  potential  of  i 

°yi  =  standard  chemical  potential  of  i 
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w-j  =  partial  molar  area  of  i  in  boundary 
ft  =  area  of  boundary 
o  =  surface  tension  of  boundary 

^  =  Helmholtz  equivalent  of  chemical  potential  of  i 
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